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FOREWORD 

The  Acoustics  and  Seismtcs  Laboratory  of  the  Institute  of  Science  and  Technology  at  The 
University  of  Michigan  has  been  conducting  research  on  seismic  wave  propagation  during  the 
past  two  years  under  Air  Force  Contract  AF  49(638)-911.  The  collections  of  a  considerable 
number  of  magnetic-tape  recordings  of  seismic  events  under  this  and  other  contracts  have 
presented  a  problem  of  data  analysts.  The  same  problem,  only  on  a  much  larger  scale,  is 
currently  being  encountered  in  the  VELA  UNIFORM  program. 

It  was  thought  that  the  use  of  analog  and  digital  computers  could  conceivably  speed  up  and 
improve  the  analyses  of  seismic  recordings.  Hence,  the  assistance  of  the  Analog  Computer 
t.aboratory  at  the  Institute  of  ,.nd  Technology  was  enlisted  to  determine  the  feasibility 

of  using  analog-computer  techniques  to  solve  many  of  the  problems. 

This  report  covers  the  theoretical  background  and  the  results  and  conclusions  that  have 
been  reached  to  date  on  the  use  of  analog-computer  techniques  for  the  analysis  of  seismic 
data  at  The  University  of  Michigan. 

The  authors  would  like  to  acknowledge  the  assistance  of  the  following  members  of  the 
Acoustics  and  Seismics  Laboratory:  J.  M.  DeNoyer,  P.  L.  Jackson,  J.  A.  Prout,  and  D.  E. 
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SYMKHS 


t 

f(t) 

Ao-  Bn 

T 

P(t) 

P 

av 

P(n) 

a(u)),  bM 

w 

W(u> 


=  Time 

=  Time  function 
«  Fourier  aeries  coefficients 

=  Fundamental  angular  frequency 


Fundamental  frequency  I— 


J' 

W^(w) 


w  .(HI,  t  ) 

u  O 


P(w.  t^) 


P(u.,to) 


(I) 


Period  of  periodic  wave,  or  length  of  portion  of  f(t)  being  analyzed 
=  Instantaneous  power  dissipation  In  1-ohm  resistor 
=  Average  power 

=  Power  In  the  n-th  trigonometric  Fourier  component 
=  Fourier  Integral  coefficients 

=  Total  energy  dissipated  across  a  1  -ohm  resistor  with  applied  voltage  f(t) 

=  Energy  contained  In  all  frequencies  from  0  to 

=  Energy-density  spectrum  (function  of  u))  (energy/(rad/sec)-bandwfdth  at  o)) 
=  Carrier  frequency 
=  Modulating  frequency 

=  Time-varying  energy-density  spectrum  at  time  I 


Double-sided  power  density  spectrum  at  time  t 


=  One-sided  power  density  spectrum  at  time  t 


2[p(u,.  t  j  = 


P(w,  t  ) 
o 


G(cu  >  (j)  ) 
o 

P(n) 

P(nuij) 

Q(t  -  t  ) 
o 

f,(t-t ) 

1  o 

P(W) 

K 

F(W) 

F(w) 


=  Filter  transfer  function 

-  Power  In  the  n-th  exponentially-defined  Fourier  component;  P(n)  =  2p{n) 

=  Power  density  In  watts  per  radian  per  second,  at  the  value  of  cu  which  Is  eq  i  -  to  nu)^ 

=  Weighting  function 

2 

=  Time  aperture  or  window:  Q(t  -  t  )  =  f,  (t  -  t  ) 

0  1  o 

°  Power  density  In  watts  per  radian  per  seeond  at  u) 

°  Normalizing  factor 
=  Fourier  transform  of  f(t) 

=  Complex  conjugate  of  F((p) 

>  Interrogation  frequency  used  in  analog-computer  circuit 
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n 

a 
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-  Response  of  an  electrical  network  to  a  unit  impulse 
•  Analog-computer  voltage 

-  Reciprocal  of  first-order  filter  time  constant 

-  Exponentially  mapped  past  (e.m.p.)  power  spectrum 


u 

o 


t 
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>  Interrogation  angular  frequency  used  in  analog-computer  circuit,  ((U^  »  2ir  F^) 

=  Frequency  coordinate  of  time-varying  power  spectrum  p(u)  ,  t  ) 

0  o 

>  Time  coordinate  of  time-varying  power  spectrum  p(ca>  ,  t  ) 

o  o 


ANALOG  COMPUTATION  OF  TIME-VARYING  POWER 
SPEaRA  OF  SEISMIC  WAVES 


ABSTRACT 

An  analog  computing  technique  for  measuring  time-varying  power  spectra  of  seismic 
signals  has  been  developed.  The  signal  to  be  analyzed  is  recorded  on  a  magnetic-tape  loop 
so  that  it  nmy  be  played  back  repeatedly,  and  the  analog  computing  equipment,  connected  to 
the  output  of  the  tape-playback  apparatus,  measures  as  a  function  of  time  the  power  level  at 
a  different  specific  frequency  for  each  passage  of  the  tape. 

The  concept  of  a  time-varying  spectrum  is  developed  and  the  theory  of  the  analysis  method 
used  is  derived.  It  Is  shown  that  in  measuring  time-varying  power  spectra,  there  is  a  theoret¬ 
ical  relationship  between  the  time  resolution  and  frequency  resolution  such  that  any  adjust¬ 
ment  which  gives  an  improvement  in  one  results  in  a  degradation  of  the  other.  This  conflict 
is  not  caused  by  equipment  limitations,  but  instead  is  a  result  of  theoretical  relations  between 
the  time  and  frequency  domains. 

A  plot  of  the  energy-  or  power-density  spectrum  of  a  seismic  wave  is  not  smooth,  but 
Instead  shows  a  number  of  sharply  defined  peaks.  It  was  found  experimentally  that  a  filter 
bandwidth  of  2  1/2%  for  the  analyzing  apparatus  gave  satisfactory  frequency  resolution  with¬ 
out  seriously  degrading  time  resolution  when  magnetic-tape  seismograms  were  analyzed. 

It  is  concluded  that  measuring  methods  based  on  the  use  of  filters  are  theoretically 
equivalent  to  methods  using  weighted  Fourier  analysis,  and  that  the  choice  of  a  method  to 
be  used  should  be  based  primarily  upon  equipment  requirements.  The  analog-computer  setup 
makes  a  very  flexible  tool  for  developing  methods  of  power-spectra  measurements,  although 
more  specialized  apparatus  would  probably  be  more  satisfactory  after  an  optimum  method 
has  been  found. 

There  Is  a  possibility  for  considerable  improvement  in  the  measuring  techniques,  and  it 
is  recommended  that  research  directed  toward  optimizing  the  analysis  of  seismic  waves  be 
continued. 


1 

INTRODUCTION 

This  report  describes  analog  computing  methods  which  were  employed  in  measuring  the  time- 
varying  frequency  spectra  of  seismic  waves.  Formal  analytical  methods,  using  the  Fourier  series 
or  Fourier  integral  formulas,  are  not  directly  applicable,  since  these  methods  yield  sinusoids  of 
constant  amplitude.  In  fact,  a  sinusoid  whose  amplitude  or  envelope  varies  as  a  function  of  time  does  not 
really  represent  a  single  frequency,  but  Instead  is  made  up  of  a  carrier  frequency  plus  sidebands.  If 
time  variations  in  amplitude  of  a  frequency  component  are  permitted,  the  value  of  frequency  cannot  be 
specified  exactly,  or  if  the  exact  frequency  of  a  component  Is  specified,  then  variations  in  its  amplitude 
are  impossible. 
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The  concept  of  a  signal  made  up  of  amplitude-modulated  frequency  components  is  not  difficult  to 
understand.  An  example  is  the  group  of  all  of  the  radio  signals  in  the  broadcast  band.  These  consist  of 
amplitude-modulated  carrier  frequencies,  spaced  10  kc  apart  and  occupying  the  frequency  band  from  550 
to  1500  kc.  A  broadcast  radio  receiver  serves  as  an  "analyzer"  in  that  it  selects  any  desired  frequency 
and  through  a  "detection"  process  produces  a  signal  corresponding  to  the  time  variation  of  the  frequency 
component's  amplitude. 

The  various  frequency  components  of  seismic  surface  waves  propagate  at  different  velocities  and,  in 
the  case  of  body  waves,  by  different  paths.  Consequently  the  arrival  timcts  for  the  different  phases  at  a 
recording  point  which  is  far  removed  from  the  source  of  the  initial  disturbance  will  occur  at  different 
times,  and  portions  of  the  wave  train  will  be  dispersed. 

In  analyzing  the  seismogram,  it  is  Important  to  determine  accurately  the  time  of  arrival  and  duration 
of  each  major  frequency  component.  Any  frequency-analysis  method  employed  should  be  capable  of  ac¬ 
curately  measuring  the  time  variations  of  a  specific  frequency  component,  but  should  at  the  same  time  be 
capable  of  separating  the  various  frequency  components.  These  two  requirements  conflict,  since  improve¬ 
ment  of  performance  on  one  of  them  results  in  degradation  In  performance  on  the  other.  The  following  • 

section  of  this  report  discusses  the  basic  problem  and  Its  solution  in  more  detail. 

2 

THEORY 

2.1.  CONVENTIONAL  SPECTRAL  ANALYSIS 

2. 1.1.  FOURIER  SERIES.  If  a  seismic  signal  is  defined  for  the  interval  0  <  t  <  T,  it  can  be  synthesized 
within  this  interval  from  a  Fourier  trigonometric  series.  The  terms  of  the  Fourier  series  are  harmonically 
related.  If  the  length  of  the  signal  being  analyzed  is  T  seconds,  the  lowest  or  fundamental  frequency  of  the 
Fourier  series  expansion  would  * j  “  j  All  other  frequencies  of  the  Fourier  terms  would  then  be 

integral  multiples  of  this  fundamental  frequency.  Thus,  within  the  interval  0  <  t  <  T,  the  signal  would  be 
represented  by  a  series  of  the  form 


cos  nu)  t  +  B 
1  n 


sin  nuijt) 


where  0)^  =  27rf ^  =  2jr/T  (radians  per  secomi)  and  the  coefficients  are  defined  as  follows: 


A 


o 


f(t)dt 


(1) 


(2) 


2 
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A 

n 


B 

n 


f(t)  cos  nu)^tdt 


f(t)  sin  nu)^tdt 


(3) 


(4) 


If  the  Fourier  series  for  a  function  f(t)  is  computed  and  the  terms  are  added  as  specified  by  Equation  1, 
the  result  is  f(t)  for  0  <  t  <  T,  but  the  function  is  periodic.  That  is,  f(nT  +  t)  =  f(t)  for  (n  =  0,  ±1,  ±2,  .  .  .). 
Thus,  the  function  obtained  from  the  Fourier  series  does  not  exactly  correspond  to  the  original  function, 
for  which  f(nT  +  t)  =  0  for  (n  =  ±1,  ±2,  ±3,  .  .  .). 

Plots  of  A  vs.  n,  and  B  vs.  n  result  in  "line"  spectra,  which  have  values  for  Integral  values  of  n  only, 
n  n 

Plots  of  this  type  are  shown  In  Figure  1. 


FIGURE  1.  EXAMPLES  OF  "LINE"  FREQUENCY  SPECTRA 


If  f(t)  is  a  voltage  which  is  applied  across  a  1-ohm  resistive  load,  the  instantaneous  power  dissipated 
in  the  load  at  any  time  t  is  given  by 

P(t)  =  f^ft)  (5) 

and  the  average  power,  in  the  interval  0  <  t  <  T,  is  given  by 

1  2 
Pav--|  If(t)l  dt 

2 

Since  the  average  power  produced  in  the  load  by  a  sinusoidal  vi  'tage  of  amplitude  A  is  A  /2.  the  average 
power  level  of  the  n-th  harmonic  of  the  Fourier  series  is  given  by 


(H) 


3 
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The  sum  of  the  average  power  levels  of  the  Fourier  components  is  equal  to  the  average  power  in  f(t) 
in  the  interval  from  0  <  t  <  T. 

Thus 


E 


2 


n=0 


|f(t))^dt 


(7) 


A  power  spectrum  for  f(t)  may  be  obtained  by  plotting  P(n)  as  a  function  of  n.  A  sketch  of  a  typical 
plot  is  shown  in  Figure  2. 


n  FREQUENCY 

FIGURE  2.  EXAMPLE  OF  A  LINE  POWER  SPECTRUM 

Note  that  the  function  P{n)  has  values  only  for  integral  values  of  n.  The  value  of  P(n)  for  each  n  cor¬ 
responds  do  the  power  level  at  that  particular  harmonic  of  the  fundamental  frequency. 

Note  also  that  the  average  power  level  P(n)  for  any  frequency  component  n  is  not  a  function  of  time. 

Thus  with  conventional  Fourier  analysis,  the  concept  of  a  frequency  component,  whose  aveiage  amplitude 
varies  as  a  function  of  time  is  not  applicable. 

2.1.2.  FOURIER  INTEGRALS  Fourier  integral  analysis  may  be  used  to  describe  the  frequency  spectra 
of  nonperiodic  waveforms.  The  Fourier  integral  formulas  are  derived  from  the  Fourier  series  formulas 
by  letting  the  period  of  the  waveform  under  analysis  approach  infinity.  This  gives  infinitesimal  frequency 
separation  between  the  frequency  components,  and  thus  leads  to  a  dense  or  continuous  frequency  spectrum, 
as  opposed  to  the  "line"  spectra  obtained  from  the  Fourier  series  expansion. 

The  Fourier  integral  representation  of  a  fun  ilon  f(t)  is  [  1)  as  follows 

r»  .00 

f(t)  -  I  a(u))  cos  uitdui  f  I  b(iii)  sin  uitdu)  (8) 

*'o  Jo 


4 


Imlituto  of  Scionco  and  Tochnology 


Tho  Univoriity  of  Michigan 


where 


a(u))  ^  I 


f(t)  cos  ultdt 


(9) 


b«iJ) 


1 

r 


f(t)  sin  uitdt 


(10) 


If  f(t)  is  a  voltage  which  is  applied  across  a  1-ohm  resistive  load,  the  total  energy  W  dissipated  for 
all  time  from  -oo  to  -mo  is  given  by 

f  *  2 

W=  lf(t)|  dt  (11) 

In  order  for  W  to  be  finite,  the  actual  duration  of  f(t)  must  also  be  finite.  If  this  requirement  is  satis¬ 
fied  (as  it  would  be  for  a  seismic  wave),  the  Fourier  Integral  energy  theorem  |1)  applies  and 


|f(t)l^dt 


■■r 

Jo 


.2 


2 


|a(ui))  |b(u>)|  dui 


(12) 


Substituting  Equation  11  into  12  gives 


W 


=  a  |a(u))| 

•'n 


|b(u))]  du> 


(13) 


Equation  13  gives  the  total  energy  in  f(t)  as  a  summation  of  the  separate  energies  of  all  the  frequency 
components.  The  energy  contained  in  all  frequency  components  between  0  and  n)  would  be  given  by 

j 


-Olj 

W(u).)  =  IT  I 

J  Jo 


(a(i»)))^  +  (b(u))|^du) 


(14) 


Since  the  integrand  is  always  positive.  W(ui^)  Is  a  monotonic  increasing  function.  This  function  may  be 
plotted  as  shown  in  Figure  3. 

Since  the  slope  of  the  curve  of  Figure  3  corresponds  to  the  Increase  In  energy  W  per  unit  Increase 
In  angular  frequency  u),  the  energy-density  function  may  be  obtained  by  differentiating  Equation  14  with 


respect  to  ui 
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FIGURE  3.  INTEGRATED  ENERGY-DENSITY  SPECTRUM 


FIGURE  4.  ENERGY-DENSITY  SPECTRUM 


'"i'  f  2  2l 

Wd(u)j)  = -^  =  ir|la(u))!  +IbM|  I  (15) 

dW(u).) 

A  plot  of  — - — -  corresponding  to  the  W(u»,)  function  of  Figure  3  Is  shown  In  Figure  4. 

d«j  J 

This  plot  corresponds  roughly  to  the  Fourier  series  power  spectrum,  as  shown  In  Figure  2,  except 
that  It  Is  a  continuous  spectrum  Instead  of  a  line  spectrum,  and  the  ordinate  Is  energy  density  Instead  of 
power. 

Thus,  the  expression 


Wdlw)  =  ir|la(u)))^  + 


(15) 


corresponds  to  energy  per  unit  bandwidth  In  the  frequency  spectrum  of  f(t).  The  energy  at  any  single 

f  2  21 

frequency  fa  zero,  since  this  Is  given  by  ir  j(a(u)))  +  Jb(w))  1  dm. 
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Aa  In  the  case  of  Fourier  series,  the  frequency  components  a(ui)  and  b(u))  are  not  functions  of  time. 

The  integral 

r® 

((t)  =  I  [a(a>)  cos  ut  +  b(iii)  sin  uit|dcc  (8) 

Jfl 

may  be  visualized  as  a  summation  of  sinusoids  of  all  frequencies,  each  sinusoid  extending  in  time  from 
-00  to  +00.  These  sinusoids  add  to  produce  f(t)  for  the  interval  in  time  during  which  f(t)  exists,  and  add  to 
produce  zero  at  all  other  times.  The  energy  in  any  one  of  these  sinusoids  is  zero,  whereas  the  energy 
density  per  unit  bandwidth  is  not,  as  shown  by  Equation  15  above.  Since  the  sinusoids  extend  in  time  from 
-00  to  -MO,  and  since  the  energy  per  unit  bandwidth  is  finite,  the  power  per  unit  bandwidth  is  zero.  If  the 
power  in  any  band  of  frequencies  were  greater  than  zero,  this  power  integrated  for  time  from  -oo  to  +oo 
would  give  infinite  energy! 

Although  Fourier  integral  analysis  can  be  applied  to  seismic  waves,  and  the  energy-density  spectrum 
can  be  plotted,  the  results  are  not  particularly  useful;  they  apply  to  the  wave  in  its  entirety  and  there  is 
no  provision  for  obtaining  the  relative  intensities  of  the  various  frequency  components  for  different  por¬ 
tions  of  the  seismogram  as  a  function  of  time.  If  we  wish  to  examine  various  parts  of  the  wave  trains  for 
their  frequency  compositions  and  plot  these  frequency  spectra  as  functions  of  time,  the  formal  techniques 
of  Fourier  integral  analysis  must  be  modified. 

2.2.  TIME-VARYING  POWER  SPECTRA 

Not  all  portions  of  a  seismogram  have  the  same  frequency  content.  That  of  the  Initial  portion  of  a  seismo¬ 
gram  may  be  completely  different  from  that  of  the  final  portion.  Thus  it  is  desirable  to  be  able  to  observe  the 
amplitude  of  each  frequency  component  as  a  function  of  time.  A  formal  Fourier  integral  analysis  of  a  seismic 
waveform  yields  a  frequency  spectrum  which  does  not  vary  with  time;  therefore,  a  modified  analytical  procedure 
is  required  to  obtain  a  time-varying  spectrum. 

2.2.1.  NATURE  OF  A  TIME-VARYING  POWER  SPECTRUM  Exact  specification  of  any  general 

function  f(t)  in  terms  of  time-varying  frequency  spectra  is  quite  difficult,  or  even  perhaps,  impossible. 

This  difficulty  results  from  the  fact  that  a  sinusoid  of  constant  frequency  whose  amplitude  is  varying  with 

time  cannot  be  regarded  as  a  single  frequency.  Instead,  since  it  is  an  amplitude -modulated  wave,  it  has 

a  frequency  spectrum  consisting  of  a  carrier  frequency  plus  bands  of  adjacent  frequencies  called  sidebands 

ll|.  If  a  sinusoid  with  a  frequency  u)  is  amplitude  modulated  by  another  sinusoid  with  a  frequency  cc  . 

c  m 

the  resulting  modulated  signal  contains  three  frequencies:  the  carrier,  ul  ;  the  upper  sideband,  ui  +  (C  ; 

c  cm 

and  the  lower  sideband,  w  -  U)  .If  the  modulation  envelope  of  the  sinusoid  is  not  itself  a  sinusoid,  then 
c  rr. 

each  Fourier  frequency  component  of  the  modulating  signal  generates  an  upper-  and  lower-sideband  fre- 
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quency.  Thus,  the  spectrum  of  the  upper  sideband  Is  the  same  as  that  of  the  modulating  signal,  with  all 
frequencies  displaced  upward  by  the  amount  of  the  carrier  frequency,  whereas  the  lower-sideband  frequency 
spectrum  is  a  mirror  image,  reflected  in  the  carrier,  of  the  upper  sideband.  These  relationships  are 
illustrated  in  Figure  5. 


(a)  (b) 

FIGURE  5.  SPECTRA  RESULTING  FROM  AMPLITUDE-MODULATION  PROCESS,  (a)  Spectra  of  modulating 
signal  and  carrier.  <b)  Spectrum  of  modulated  sinusoid. 


The  bandwidth  or  frequency  spread  occupied  by  the  signal  is  determined  by  the  highest  frequency  com¬ 
ponent  of  the  modulation  envelope.  Consequently,  if  we  expect  to  detect  rapid  variations  in  the  amplitude 
of  a  given  frequency  component  of  a  seismic  signal,  we  must  really  examine  the  signal  over  a  wide  band 
of  frequencies;  therefore  it  becomes  impossible  to  specify  exactly  what  frequency  component  is  being 
measured.  On  the  other  hand,  if  a  very  narrow  band  of  frequencies  is  being  examined,  so  that  the  fre¬ 
quency  component  whose  amplitude  is  being  measured  is  specified  very  accurately,  rapid  variations  in 
the  amplitude  of  this  component  cannot  be  detected  because  the  sideband  frequencies  which  are  generated 
by  the  higher  frequency  components  of  the  mjdulation  envelope  are  being  suppressed. 


2.2.2.  TECHNIQUES  FOR  COMPUTING  TIME- VARYING  POWER  SPECTRA 


2.2.2. 1.  The  Rectangular  Time  Aperture.  One  possible  method  of  computing  the  time-varying  power 
spectrum  of  a  signal  f(t)  is  illustrated  by  Figure  (i. 

Suppose  the  actual  duration  of  f(t)  is  unknown,  and  we  wish  to  estimate  its  frequency  composition  at 

time  t  .  The  portion  of  f(t)  extending  from  t  -  T/2  to  t  +  T/2  could  be  analyzed  to  determine  the  energy- 
o  o  0 

density  spectrum  for  this  portion  of  f(t).  The  energy-density  spectrum  for  all  of  f(t)  cannot  be  computed 
from  this  sample,  but  the  energy  contributed  during  this  period  can  be  computed  and  divided  by  the  length 
of  the  period  to  give  an  average  power-density  spectrum  for  the  Interval  from  t  -  T/2  to  t  +  T/2. 
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FIGURE  6.  EXAMINATION  OF  f(t)  THROUGH  A  "TIME  WINDOW" 


Thus,  if  we  consider  f(t)  as  existing  only  in  the  interval  (t  -  T/2)  <  t  <  (t  +  T/2),  we  can  compute 

o  o 

the  energy-density  spectrum  as  follows:  changing  the  limits  for  the  integrals  of  Equations  8  and  9,  and 

considering  the  results  as  a  function  of  t  we  derive 

o 


t  +T/2 
o 


1  f 

a(u>.  tJ  =  7 

A  -T/2 


f(t)  cos  uitdt 


(16) 


b(u»,  t 


o  X  J 


t  +T/2 


f(t)  sin  uitdt 


t  -T/2 
o 


(17) 


The  energy-density  spectrum  is  given  by  Equation  15; 


W  (u»  =  x{(a{u»|^  +  (b(u))|^} 


(15) 


and  the  average  power-density  spectrum  at  t^  may  be  obtained  by  dividing  by  the  length  of  the  interval; 

Substituting  Equations  16  and  17  into  Equation  18  gives 


P(W,  t  1  =  — 
o  xT 


f 

r  t  *T/2 

r  ° 

2 

r  t  +T/2 
fo 

1 

) 

1  f(t)  cos  uXdt 

1  + 

J  f(t)  sin  uJtdt 

1 

•'t  -T/2 

U  o  —i 

t  -T/2 
^  0 

J 

(19) 


The  integrals  in  Equation  19  could  be  evaluated  continuously  on  analog  computing  equipment  for  any 

specified  u)  as  a  function  of  t  by  using  tape-reconiing  equipment  with  spaced  playback  heads,  so  that  the 
o 
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T  T 

Integrals  of  f(t)  could  be  obtained  simultaneously  between  the  limits  of  0  and  ^ ^  ~  ~2  ' 

The  first  integral  of  Equation  19  can  be  rewritten 


t  +T/2 


t  +T/2 


f" 

I  f(t)  cos  lUtdt  =  I  f(t)  cos  uitdt  - 

•'t  -T/2  •'o  Jo 

o 


t  -T/2 


f(t)  cos  UHdt 


(20) 


providing  f(t)  is  continuous  from  0  to  t^  +  T/2.  and  t^  -  T/2  >  0. 

Substituting  t  +  T/2  for  t  in  the  first  Integral  on  the  right  side  of  Equation  20  gives 


t  +T/2  ( 

f  ° 

I  f(t)  cos  Uitdt  =1  f(t  +  T/2)  cos  u)(t  +  T/2)dt 

■'-T/2 


(21) 


Similarly,  the  second  integral  may  be  written 
t  -T/2 


J  f(t)  cos  Uitdt  =  J 


f(t)  cos  Uitdt  =  j  f(t  -  T/2)  cos  ui(t  -  T/2)dt 

■'t/2 


(22) 


If  f(t)  =  0  for  t  <  0,  then  f(t  -  T/2)  =  0  for  t  <  T/2,  and  the  lower  limit  of  the  integral  on  the  right  in 
Equation  22  may  be  extended  to  -T/2  without  affecting  the  value  of  the  result.  Therefore,  substituting 
Equations  21  and  22  in  20  gives 


/ 


t  +T/2 
0 


f(t)  cos  Uitdt 


t  -T/2 
0 


lf(t  +  T/2)  cos  ui(t  +  T/2)  -  f(t  -  T/2)  COS  ui(t  -  T/2))dt 


(23) 


Equation  23  could  easily  be  mechanized  on  analog  computing  equipment  in  combination  with  a  tape 
recorder  having  two  playback  heads  spaced  T  seconds  apart,  so  that  f(t  +  T/2)  and  f(t  -  T/2)  could  be 
obtained. 

A  similar  method  could  be  used  to  obtain 


/. 


t  +T/2 
o 


f(t)  sin  uitdt 


t  -T/2 
o 


and  F(ui,  t  )  as  a  continuous  function  of  t  could  be  computed  from  Equation  19. 
o  o 


(24) 
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» 


There  would  be  no  need  for  a  tape  recorder  if  a  device  having  a  transport  delay  of  T  seconds  were 
available.  A  block  dlagrann  showing  the  method  of  computing  the  first  integral  of  Equation  19  be  using 
such  a  device  is  shown  in  Figure  7. 


FIGURE  7.  ANALOG-COMPUTER  MECHANIZATION  OF  EQUATION  23 


2. 2. 2. 2.  The  Rectangular  Frequency  Aperture.  The  power  spectrum  of  f(t)  as  a  function  of  time  might 
also  be  found  by  use  of  a  filter.  For  example,  a  filter  having  a  finite  passband  with  sharp  cutoff  at  the 
edges  of  the  passband  might  be  regarded  as  a  rectangular  aperture  or  "window"  for  viewing  the  frequency 
spectrum  of  f(t).  Figure  8  shows  the  passband  of  such  a  filter  superimposed  upon  P(u)),  the  power  spectrum 
of  t(t). 


The  filter,  whose  pass  charcteristic  is  shown  in  Figure  8,  passes  all  frequencies  from 
Au) 

u)^  +—  without  attenuation,  and  completely  suppresses  all  frequencies  outside  this  range. 

function  is  G((i)  -  UJ  ),  the  power  transmitted  through  the  filter  is  given  by 
o 


Aw , 

w  -  to 
O  2 

If  the  filter 


p(Aw)  =  I  p(u))G(w  -  u>  )dw  = 
•^0  ° 


w 


U)  - 

o 


Aw 

'  2 

Aw 

2 


p(w)  dw 


(25) 


FIGURE  8.  ILLUSTRATION  OF  A  FREQUENCY  APERTURE 
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and  the  average  power  per  unit  bandwidth  for  the  frequency  range  covered  by  the  filter  la  obtained  by 
dividing  Equation  25  by  diO). 

Aw 
+  — 

1  c  “  2 

P(w  .  ‘)  =  ^  PM  dw  (26) 

J  Aw 

o  2 

The  output  of  the  filter  may  be  time-varying,  subject  to  restrictions  Imposed  by  the  filter's  bandpass 

characteristics.  The  filter  should  have  a  narrow  passband  In  order  to  restrict  the  range  of  w  over  which 

the  power  Is  averaged  for  any  given  w  ,  and  thus  give  good  frequency  resolution.  However,  If  the  pass- 

o 

band  Is  very  narrow,  rapid  variations  of  the  power  output  cannot  be  transmitted  because  the  sideband  fre¬ 
quencies  which  must  be  part  of  such  a  signal  are  suppressed. 

In  practice,  the  filter  would  be  set  for  various  values  of  w^,  and  the  output  power  would  be  Interpreted 

as  Aw  P(W  ,  t). 
o 

2. 2, 2. 3.  Comparison  of  Time-Aperture  and  Frequency- Aperture  Approaches.  Summarizing  Sections 
2.2.1  and  2.2.2,  It  appears  that  there  are  two  possible  methods  of  obtaining  time-varying  power-density 
spectra  of  f(t).  One  of  these  Involves  examination  of  a  portion  of  f(t)  through  a  moving  "window"  In  time 
and  performing  Fourier  Integral  analysis  on  the  portion  of  f(t)  seen  through  the  window.  The  second  In¬ 
volves  use  of  a  filter,  or  frequency-spectrum  window,  through  which  the  power  content  of  a  band  of  fre¬ 
quency  components  of  f(t)  are  observed.  The  entire  frequency  spectrum  could  be  covered  either  by  using 
a  number  of  filters,  or  by  adjusting  the  center  frequency  of  a  single  filter  and  making  consecutive  obser¬ 
vations. 

Increasing  the  length  of  the  time  window  for  the  first  method  improves  the  accuracy  of  frequency- 
component  determination,  but  reduces  the  rate  at  which  the  power  density  being  measured  can  vary.  In¬ 
creasing  the  length  of  the  frequency  window  for  the  second  method  has  the  opposite  effect.  It  increases 
the  rate  at  which  the  output  of  the  filter  may  vary  but  reduces  the  selectivity,  or  frequency  resolution. 

It  will  be  shown  later  that  the  two  methods  discussed  above  are  not  really  separate  and  that  any  method 
of  analyzing  a  signal  f(t)  In  terms  of  a  time-varying  frequency  spectrum  involves  the  use  of  both  a  time 
aperture  and  a  frequency  aperture. 

2.2.2. 1.  The  Variable-Transmission  Time  Aperture.  Use  of  the  rectangular  time  aperture,  as  des¬ 
cribed  in  Section  2.2.1,  gives  the  average  power-density  spectrum  of  f(t)  for  the  interval  (t  -  T/2)  <  t  < 

o  “  ~ 

(t  +  T/2).  In  determining  this  average  P(w,  t  ),  equal  weight  Is  given  to  the  spectral  characteristics  of 
o  o 

f  (t)  throughout  the  Interval  covered  by  the  time  aperture.  Better  time  resolution  would  result  If  the  spec- 
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tral  characteristics  of  f(t)  in  the  neighborhood  of  t^  were  emphasized  in  the  computing  process.  This  may 

be  accomplished  by  using  a  variable-transmission  time  aperture  whose  transmission  is  maximum  in  the 

neighborhood  of  t  .  An  equation  for  conq>uting  P(ii),  t  ),  which  uses  a  time  aperture  of  this  type,  will  be 
o  o 

derived  in  this  section. 

In  deriving  this  formula,  it  is  comvenient  to  use  the  exponential  form  of  the  Fourier  series  and  inte¬ 
gral  formulas.  For  the  Fourier  series,  these  are  (21 


oo 

f(t)  =  F(n)e 
n=-oo 


jnWjt 


(27) 


F(n) 


1 

=  T  f(t)e 

J-T/2 


-jnWjt 


dt 


(28) 


F(n)  =  ("  =  *1'  *2 - )  (29) 

4  n  n 

where  A  and  B  correspond  to  the  coefficients  of  the  trigonometric  series  (Equation  1).  The  factor  of 
n  n 

1/2  occurs  in  Equation  29  because  the  exponential  form  of  the  series  uses  both  negative  and  positive  values 
of  n,  and  thus  has  twice  as  many  terms  as  the  trigonometric  series,  which  uses  only  positive  values  of  n. 

If  f(t)  is  the  voltage  across,  or  the  current  through  a  1  -ohm  resistor,  the  average  power  dissipated  in 
the  resistor  by  the  n-th  harmonic  component  is 


.2  „  2 
A  +  B 


p(n)  =  lF(n)|  =  — - 2-  (n  -  o,  ±1,  ±2,  .  .  ,) 

•f 


(30) 


Comparison  of  Equation  30  with  Equation  6  reveals  that  the  p(n)  of  the  former  is  equal  to  half  the  p(n) 
of  the  latter.  Similarly,  any  time-varying  continuous  power  spectrum  derived  from  the  exponential  for¬ 
mulas  will  include  both  positive  and  negative  values  of  u),  and  therefore  the  value  of  p(u)  ,  t  )  so  obtained 

o  o 

must  be  doubled  to  obtain  P(u)  ,  t  ),  the  physical  power  level  at  frequency  IP  at  time  t  . 

0  0  0  0 

Substituting  Equation  28  into  Equation  30  gives 
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where 

„  ,  2ir 

(A  =  2irf.  =  — 

1  1  T 

and  is  the  fundamental  frequency. 

Equation  31  gives  the  power  level  of  each  frequency  component.  Since  each  harmonic  frequency  Is  an 
Integral  multiple  of  the  fundamental  frequency,  a  plot  of  p(n)  versus  n  would  be  a  "line"  spectrum,  with  Wj 
as  the  spacing  between  adjacent  lines. 

Even  If  f(t).  Is  nonperlod'c,  the  Fourier  components  for  a  section  of  f(t)  extending  from  -T/2  to  T/2 
are  still  given  by  Equation  31.  The  length  of  the  section  under  consideration  may  be  extended  by  Increasing 
T.  When  T  Is  Increased,  Is  decreased,  thus  reducing  the  spacing  of  lines  on  the  power  spectrum.  As  T 
Is  allowed  to  approach  Infinity,  approaches  zero  and  the  power  spectrum  becomes  a  continuous  plot.  Since 
there  are  now  an  Infinite  number  of  frequency  components  In  any  finite  bandwidth,  and  the  power  In  such  a 
bandwidth  Is  also  finite,  the  power  per  frequency  component  must  be  zero. 

Thus. 


11m 

T-oe 


p(n)  = 


11m  1_ 

T-oo_2 

T 


(32) 


Thus,  the  power  spectrum  can  no  longer  be  expressed  In  terms  of  power  per  frequency  component.  In 
order  to  avoid  this  difficulty,  the  power  density  at  any  given  frequency  may  be  expressed  In  terms  of  power 
per  unit  bandwidth. 

If  p(n)  =  power  per  frequency  component,  and  the  frequency  components  are  spaced  f^  cps  apart,  the 
number  of  frequency  components  per  cps  Is  l/fj’  ^^d  the  number  of  frequency  components  per  radian  per 
second  Is  1/w  =  l/2irf  .  The  power  density  In  watts  per  radian  per  second,  then.  Is  p(nu),)  . 

11  1  2irlj 

Substituting  this  Into  Equation  31  gives 


or,  since  f ^  =  1/T, 


1  -J"“l‘  2 

P(nw )  - - -  f(t)e  dt 

2irfjT^  I  •^-T/2 


p(nWj)  = 


1 

2x7 


-JnW  t 

f(t)e  dt 


2 


(33) 
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Then,  in  order  to  obtain  a  continuous  frequency  spectrum,  let  T  become  infinite.  The  quantity  n(0^, 
becomes  ut  and 

IT 

f  f(t)e-l“‘dt  (34) 

J.T/2 

At  this  point,  a  brief  discussion  of  the  nature  of  the  function  defined  by  Equation  34  is  in  order.  It  may 
be  shown  that  this  function  is  not  valid  as  a  statistical  estimate  of  the  power  spectrum  of  a  random  process, 
because  the  variance  of  the  p(ia)  defined  in  this  manner  does  not  approach  zero  as  T  becomes  infinite,  ( 3] . 
Nevertheless,  Equation  34  accurately  defines  the  power  spectrum  of  the  particular  function,  f(t).  to  which 
the  equation  is  applied  in  that  the  power  spectrum  so  defined  corresponds  exactly  to  the  spectral  density  of 
the  Fourier  series  for  f(t).  The  function,  p((ii),  as  defined  by  Equation  34  is  extremely  irregular.  Sincep(U>) 
is  proportional  to  the  amplitudes  of  the  Fourier  series  components  making  up  f(t),  and  since  f(t)  is  infinitely 
long,  the  increment  between  terms  of  the  series  is  infinitesimal,  and  thus  any  small  interval  of  u)  (no  matter 
how  small)  contains  an  infinite  number  of  different  values  for  p(u>).  If  p(u>)  is  averaged  over  a  small  incre¬ 
ment  of  01,  however,  a  consistent  usable  value  for  average  p(w)  for  that  interval  is  obtained.  According  to 
Goldi  ran,  (4|,  this  averaged  or  smoothed  power  spectrum 

0>+Au) 

'  A«"lo  h>  /  ‘‘‘P 

01 

corresponds  to  the  statistical  power  spectrum  for  the  function  f(t).  In  other  words.  p(oi)  as  defined  above  is 
also  the  result  of  averaging  p(u>)  from  Equation  34  for  an  infinite  number  of  random  functions  of  time,  each 
having  the  same  statistical  power  spectrum.  In  the  steps  which  follow,  since  the  sections  of  f(t)  being  exam¬ 
ined  are  of  finite  length,  an  averaging  with  respect  to  oialways  occurs,  and  we  are  not  bothered  by  a  vari¬ 
ance  of  p(w)  of  the  type  discussed  above. 

Next,  let  us  consider  the  problem  of  estimating  p(u))  by  observing  a  finite  section  of  f(t).  If  f(t)  has  uni¬ 
form  power  spectral  characteristics  throughout  the  length  T  of  the  section  being  observed. 


2 

(35) 

for  <4)  >  >  — 

T 

Since  T  is  now  finite,  we  are  really  formally  justified  in  computing  the  frequency  spectrum  only  at 
values  of  ui  which  are  multiples  of  the  fundamental  frequency  =  2ir/T.  Evidently,  the  resolution  of  the 
frequency  spectrum  obtained  from  Equation  35  is  of  the  order  of  1/T  cps.  Roughly,  about  1%  resolution 
could  be  obtained  if  T  corresponded  to  100  cycles  of  a),  the  frequency  for  which  the  power  density  is  to  be 


p(u») 


1 

2irT 


.T/2 


f(t)e'^‘*^dt 


'-T/2 
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computed.  If  p(b))  Is  to  be  determined  accurately  at  low  frequencies,  long  Intervals  of  f(t)  must  be  observed; 
but  for  higher  frequencies,  the  interval  T  can  be  shortened. 

If  a  portion  of  the  time  axis  of  f(t)  is  divided  into  a  number  of  intervals,  each  long  enough  to  permit  a 
reasonably  accurate  computation  of  p(a>)  using  Equation  35,  and  if  the  power  density  spectrum  of  f(t)  is  uni¬ 
form  throughout  the  portion  chosen,  then  the  values  of  p(a>)  obtained  by  applying  Equation  35  to  each  of  the 
intervals  should  be  approximately  the  same. 

For  each  interval,  extending  from  time  t  ,  to  t  ,  Equation  35  may  be  written 

n-1  n 


or 


p(u>)  = 


2ir(t 

n 


) 


t 


2 


(36) 


2ir(t  -  t  )  p(a))  = 
n  n-l 


t 


2 


(37) 


Summing  all  equations  of  this  type  from  n  =  1  to  n  -  N,  where  °  assuming  that  p(u>)  is 

constant  over  this  Interval,  we  derive 

N  N 

n=I  n-l 


i 

I 


f(t)e’^‘^dt 


n-l 


(38) 


However, 

£“n-‘n-l>  =  '' 
n=l 


(39) 


therefore. 


2vp{(ji)  T 


f(t)e‘^‘^dt 


2 


(40) 


And  substituting  Equation  40  into  Equation  35  gives 


16 


Imtitut*  of  Scionco  and  Tochnelofly 


Tho  Univortity  of  Michigan 


An  integral  over  any  Interval  may  be  expresaed  as  the  sum  of  integrals  over  all  subintervals  contained  in 
this  interval;  thus 


The  operation  of  squaring  the  absolute  values  of  both  sides  of  Equation  42  and  substituting  the  result 
into  Equation  41  gives  the  rather  interesting  relation 


N 


E 

n=sl 


I 


f(t)e'^‘^dt 


n-1 


(43) 


If  (t  -  t  ,)  =  At  (with  all  intervals  equal),  and  thus 
n  n-1 

N 

V  (t  -  t  ,)  «  T  =  NAt 

Z_.  " 

n=l 


then  Equation  38  can  be  rewritten  to  give 


P(«>)  = 


E 

n=  1 


1 

2irAt 


r  f(t)e-J‘^dt 
•'t  , 


(44) 


Thus,  this  gives  p(u>)  for  the  interval  from  t  =  0  to  t  =  NAt  as  the  average  of  the  values  of  p(u))  computed 

for  the  subintervals  t  ,  <  t  <  t  . 

n-1  n 

If  the  power  spectrum  of  f(t)  varies  with  time.  Equation  44  specifies  the  average  value  of  p(u))  for  the 

section  of  f(t)  being  examined.  In  order  to  define  p(u>)  at  some  particular  time  t^,  a  weighting  function  may 

be  used  in  computing  this  average.  Suppose  we  use  a  weighting  function  Q(t  -  t  ),  which  reaches  a  maximum 

0 

for  some  particular  instant  t  and  approaches  zero  for  values  of  t  which  are  significantly  removed  from  t  , 

00  ° 
in  such  a  manner  that  the  integral  /  Q(t  -  t  )dt  exists.  (Figure  9  shows  a  typical  weighting  function.)  Use 

-oo  ° 

of  such  a  weighting  function  will  give  more  weight  to  the  spectral  characteristics  of  f(t)  in  the  neighborhood 

of  t  =  t^,  and  thus  the  weighted  average  of  p(u)),  obtained  by  using  this  weighting  function,  can  be  interpreted 

as  p(ui)  at  t  =  t  . 

o 
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FIGURE  9.  A  WEIGHTING  FUNCTION 


Using  the  weighting  function  Q(t  -  t  ),  Equation  44  can  be  rewritten 

o 


p(a),  y  = 


n-1  n-1 


N 

z 

nsl 


where  Q(t  -  t  )  la  the  average  value  of  the  weighting  function  in  the  interval  t  <  t  <  t  . 
no  n*l  n 


Since  Q(t  *  t  )  is  constant  for  any  particular  t  ,  Equation  45  can  be  rewritten 
no  n 


(45) 


P(w.  t  )  = 
o 


n-1  I  n-1 


^Q(tn-to)At 

n=  1 


(46) 


If  N  is  permitted  to  increase  now,  so  that  Equation  46  includes  all  time  from  -ao  to  oo.  Equation  46  can  be 
rewritten 


p(u),  t^)  = 


*  1 


i 


f(t)e '  dt 


2]  ‘*'‘n-‘o>^‘ 


(47) 
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and,  when  Equation  -11  la  used. 


p((iJ.  t  )  = 
o 


C>^ 


t  )  f(t)e‘^‘^dt 
0 


00 


i; 


n=-o0 

Finally,  If  we  allow  At  In  the  denominator  of  Equation  48  to  approach  zero, 

|2 


PiW.  t^)  = 


1  , _ 

if 

*'“00 


y  f(t)e'^‘*^dt 


OO 

/  «-•, 


-  t  )dt 
o 


and,  If  we  redefine  yQ(t  -  t^)  =  f^ft  -  t^), 


p(u>,  1^)  = 


1  r* 

r.  V*- 

«^“0O _ 


t  )f(t)e'^‘**dt 
0 


f 

-to 


-  t  )dt 
o 


(48) 


(49) 


(50) 


For  any  specified  weighting  function,  the  denominator  of  the  right  side  of  Equation  50  Is  a  constant: 

to 


I  f>-Vdt=i 


(51) 


So,  this  equation  can  be  rewritten 


/x> 

fj(t  -  t^)((t)e-^‘^dt 
TO 


(52) 


where 


/  C'.-, 


-  t  )dt 
o 


(53) 


Equation  52  can  be  Interpreted  as  describing  the  use  of  an  aperture  or  window  In  time,  having  a  trans- 
slon  funct 
at  that  point. 


mission  function  of  f,  (t  -  t  )  to  examine  f(t)  about  the  Instant  t  and  to  determine  the  power  spectrum  of  ((t) 
1  o  o 
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If  f,  (t  -  t  )  s  1  for  all  values  of  t  and  t  ,  Equation  52  becomes  equivalent  to  Equation  34,  which  is  restated 
1  o  o 


here: 


p(w)  = 


iim 
T-^  2tT 


r' 

^  J-T/2 


t(t)e"^‘^dt 


(34) 


Thus  we  obtain  the  non-time-varying  power-density  spectrum,  where  p(a))  is  the  average  power  density  at  w 
for  all  time,  from  -oo  to  +<o. 

Evidently,  the  time  resolution  of  p(a),  t  )  as  computed  from  Equation  52  is  improved  as  the  interval 

o 

covered  by  the  window  function  f,  (t  -  t  )  becomes  narrower,  but  a  resonably  long  window  function  must  be 

1  o 

used  to  give  good  frequency  resolution. 

2. 2.2.5.  The  Frequency-Aperture  Equivalent  of  a  Time  Aperture.  Next,  we  show  that  use  of  a  time  aper¬ 
ture  in  the  manner  specified  by  Equation  52  is  actually  equivalent  to  using  a  window  in  the  frequency  domain. 

For  power  at  a  specific  time  t  and  at  a  specific  frequency  u)  ,  Equation  52  is  written 
o  0 


P<W  .  t  ) 
o  o 


■^\I> 


-jU)  t 

t  )  f(t)e  °  dt 
o 


(54) 


Substituting  T  +  t  for  t  in  the  integrand  of  Equation  54  gives 
0 


p<w  » t  )  •  — 
'  o  o  2ir 


t 

0  0  ^ 


L 


-ju>  T 

f(T  +  t  )  f, (T)e  °  dr 
o  1 


(55) 


Since  (0^  is  a  constant  for  this  integration,  we  can  define 


fgfT)  '  fj(T)e 


-Jut  T 
O 


(56) 


Substituting  this  in  Equation  55  gives 


p{U)  ,  t  )  “  ~  p 
o  o  2ir  f 


-JU>  t  .03 


J-OO 


(57) 


However,  it  can  be  shown  ( 2)  that 


]iut 

f  (T)f(T  + 1  )dT  »  I  27r  F  ((i;)F(w)e  °du) 


f  yT)f(T*t  )dT-  r 

J-  QO  J-OB 


(68) 


where 


F(W) 


2ir  I 

%/-C0 


f(T)e  ^**^dT 


(69) 
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Equation  61  can  be  rewritten 


Vr)e^‘^dr 

J-QO 

J  -aa 


(60) 


(61) 


Since  by  definition. 


-  fj(T)e 


JU)  T 

o 


(56) 


1 

=2tJ  S<- 


^oo  j(aHii)  )t 

F2(w)=srl  f,(^)e  °  dT 


=  F,(W  +  ti> ) 
1  0 


(62) 


Substituting  Equation  (>2  into  Equation  58  and  the  result  into  Equation  57. 


p(«  .  t  )  =  _ 
0  o  2v 


-III)  t  r* 


—  iu,*t 

2)rF(ii))F,(ii)  +  l»)  )e  du) 
1  o 


(63) 


-Jo)  t 
o  o 

and,  since  e  -  1 


p(<D__ 


t  )  =  2irK 
o 


If 

M-30 


—  Ju^ 

F(u))F,(u>  +  u>  )e  du> 
1  o 


(64) 


The  factor  K  In  Equation  64  can  be  rewritten  In  terms  of  F,((i)  +  u)  )  by  using  the  Fourier  integral 

1  o 

energy  theorem. 

Since  f,  (t)  and  F,  (u)  +  u)  )  are  a  Fourier  transform  pair, 

1  1  o 


_1_ 

K 


•'-00  •'-OO  ■'-00 


2ir 


Fill)  +  (i) ) 
1  o 


dll) 


(65) 


and  since 


F,(u)  +  ii))  =  F,(li)  +  li)) 

1  O  10 


(('>(>) 
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therefore 


K  = 


2rJlFj(U.  +  a. 


)  ^da> 


m 


Substituting  Equation  <>7  into  Equation  (>4.  substituting  -u)  for  ti)  in  the  result,  and  recalllngthat  p((i>  ,  t  )  = 

O  O  0  0 

p(-(j  I  t^)  because  of  symmetry  In  the  power  spectrum,  we  obtain; 


0  0  o  o 


r“ 

L 

Fjfll 

r 

L 

I )  F(u))e  '’diu 


(fiS) 


dill 


Repeating  Equation  SO  for  comparison. 


P(<‘>„.  t  )  = 
o  o 


-Ju)  t 

f,(t-t)f(t)e  °  dt 
1  o 


2 


t  )dt 
o 


(50) 


In  comparing  Equations  50  and  <;8.  we  see  clearly  that  since  Equation  50  describes  p(a>  ,  t  )  In  terms 

0  0 

of  a  portion  of  fit)  as  seen  through  a  variable-transmission  "time-window”  function  fj(t  -  t^).  Equation  68 

describes  P((i>  ,  t  )  in  terms  of  a  portion  of  F(a)),  the  spectrum  of  F(t)  as  seen  through  a  varlable-trans- 
o  o  __  _ 

mission  "frequency-window"  function  F,(<ii  -  a> ).  The  frequency  window  F,(ii»  -  U) )  is  the  complex  con- 

10  1  o 

jugate  of  the  Fourier  transform  of  the  time  window,  f  (t  -  t  ). 

1  o 

In  general,  if  the  time  aperture  is  narrowed  in  order  to  obtain  better  time  resolution,  the  effective 
filter  bandwidth  is  broadened,  and  the  frequency  resolution  becomes  poorer.  If  the  time  aperture  is  widened 
so  that  the  frequency  spectrum  is  computed  for  larger  sections  of  f(t),  the  frequency  resolution  is  improved, 
but  the  time  response  or  resolution  is  degraded.  Specification  of  either  the  time  aperture  or  the  filter  band¬ 
width  characteristic  determines  both,  since  one  is  derived  from  the  other.  The  exact  choice  of  either 
aperture  depends  upon  the  nature  of  the  data  to  be  analyzed  and  the  type  of  spectral  information  desired. 

2.2.3.  THE  FREQUENCY-ANALYSIS  METHOD  USED  IN  THE  1ST  ANALOG  COMPUTER  LABORATORY 
2,2.3. 1.  Computer  Setup.  In  preparation  for  measurement  of  its  power  spectrum  on  the  analog  com¬ 
puter,  the  seismic  signal  to  be  analyzed  was  first  recorded  on  one  channel  of  a  length  of  magnetic  tape. 

The  ends  of  this  length  of  tape  were  spliced  together  to  form  a  loop  so  that  when  the  loop  was  played  back, 
the  recorded  seismic  signal  would  be  reproduced  repeatedly.  Control  signals  were  recorded  on  a  second 
channel  of  the  tape  at  positions  corresponding  to  the  beginning  and  ending  of  the  seismic  signal.  These 
control  signals  were  used  to  program  the  operation  of  the  analog  computer. 
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The  computer  was  so  mechanized  that  It  examined  a  different  segment  of  the  power  spectrum  of  the 
seismic  signal  for  each  pass  of  the  tape,  and  plotted  the  results  automatically.  The  control  signals  which 
were  recorded  on  the  second  channel  of  the  tape  actuated  the  analog  computer  so  that  it  automatically  reset 
itself  to  a  new  "interrogation"  frequency  for  each  pass  of  the  tape.  Thus,  each  pass  of  the  tape  generated 
Information  concerning  a  new  portion  of  the  power  spectrum. 

Figure  10  is  a  block  diagram  of  the  basic  computing  circuit  used.  The  oscillator,  operating  at  some 
interrogation  frequency  generates  signals  corresponding  to  cos  U^t  and  sin  bi^t.  Multiplier  A  multiplies 
the  incoming  seismic  signal  f(t)  by  cos  u^t,  whereas  multiplier  B  multiplies  f(t)  by  sin  u>^t. 

The  output  signal  of  multiplier  A  (designated  on  the  diagram)  next  passes  through  filter  A,  and  the 

filtered  signal  E„  is  squared  in  multiplier  C  to  produce  E.,  which  is  fed  into  the  adder. 

3  0 

The  output  of  multiplier  B  is  similarly  processed  to  produce  a  second  Input  to  the  adder. 

2.2. 3.2.  Theory  of  Computer  Circuit  Operation.  The  operation  of  the  circuit  can  be  explained  qualita¬ 
tively  in  terms  of  heterodyning  and  filtering.  Multiplication  of  f(t)  by  cos  lU^t  shifts  its  frequency  spectrum. 
F(<i)),  by  Wj  so  that  the  frequency  spectrum  of  Ej  is  F(u)  -  Wj).  Therefore,  examination  of  the  frequency 
spectrum  of  E^  for  zero  frequency  or  d-c  components  is  equivalent  to  examination  of  f(t)  for  uij  components. 


Multiplier  A  Filter  A  Multiplier  C 


Multiplier  B  Filter  B  Multiplier  D 

FIGURE  10.  FUNCTIONAL  BLOCK  DIAGRAM  OF  ANALOG-COMPUTER  CIRCUIT  USED  IN  FREQUENCY 

ANALYSIS 
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Filter  A  Is  a  low-pass  filter  which  passes  d-c  perfectly  but  attenuates  a-c  signals.  Use  of  this  multiplier 
and  lowpass  filter  la  roughly  equivalent  to  using  a  passband  filter  centered  at  ai|  on  f(t)  directly.  The  output 
of  filter  A  Is  squared  by  multiplier  C,  thus  converting  voltage  amplitude  to  power. 

The  lower  channel  operates  similarly,  except  that  It  multiplies  f(t)  by  sin  u^t  and  produces  a  signal 
(E^)  having  a  frequency  spectrum  which  Is  in  phase  quadrature  with  that  of  (E^).  After  this  signal  is  passed 
through  lowpass  filter  B  and  squared,  it  la  added  to  the  squared  signal  from  the  upper  channel  and  the  sum 
is  then  a  voltage  which  Is  proportional  to  pfai^,  t). 

A  more  quantitative  explanation  of  the  operation  of  the  analog-con>puter  circuit  can  be  made  by  showing 
that  It  evaluates  p((U^,  t^)  In  terms  of  Equation  52,  which  was  derived  above.  Rewriting  Equation  52, 


K  r® 

P<‘"i-V'27  I 

J-aa 


(52) 


where 


A  QO 

j  fj®(t-t)dt  (53) 

J.oo 


and  where  fj(t  -  t^)  is  the  time-aperture  function 
f(t)  is  the  function  being  analyzed 

and  t^  are  the  frequency  and  time  coordinates  at  which  p  is  being  evaluated 


-JWjt 

Since  e  >  cos  u^t  -  ]  sin  (U^t,  Equation  52  can  be  rewritten 


J-OO 


t  )f(t)  COB  w.tdt  ■ 
o  1 


t  )f(t)  sin  w,tdt 
o  i 


(69) 


And,  since  the  square  of  the  magnitude  of  a  complex  quantity  is  equal  to  the  sum  of  the  squares  of  its  real 
and  imaginary  components. 


p<"i-  V 


f 

2 

f* 

2 

i  J  ‘^i‘  *** 

+ 

I  fj(t  -  yf(t)  sin  taijtdt 

J  -CO 

(70) 
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then 


<79) 


where  p(w,,  t  )  is  the  magnitude  of  the  time-varying  power  spectrum  as  defined  by  Equation  70  for  time  t 
10  o 

and  angular  frequency 


and 


/ 


oo 

oO 


1 

f,^t-t  )dt 
1  o 


In  this  last  expression 

fi<t-V-h(-(t-yi 

where  h(t  -  t  )  la  the  response  of  the  Alters  to  a  unit  Impulse  occurring  at  time  t  , 
o  o 

Thus,  the  analog-computer  circuit  evaluates  Pfo)  .  t  |  In  accordance  with  Equation  70.  and  uses  a  time 

i  o 

aperture  or  weighting  function  which  is  equivalent  to  the  fitter's  input  response,  folded  about  the  ordinate 
t  =  t^.  The  relationship  between  h(t  -  t^)  and  fj(t  -  t^)  is  shown  in  Figure  11. 

2.2.3.3.  Computer  Circuit  Operation  Using  First-Order  Lowpass  Filters.  In  Equation  79  above,  it  is 

assumed  that  the  filter  impulse  response  is  such  that  the  time-window  function,  f,(t  -  t  )  =  hi  -(t  -  t  )) , 

10  o 

gives  more  weight  to  the  characteristics  of  f(t)  in  the  neighborhood  of  t  and  less  weight  to  the  parts  of  f(t) 

o 

which  precede  t^.  Next,  we  shall  examine  the  effects  of  using  various  filter  transfer  functions. 

A  first-order  low-pass  filter  is  easily  synthesized  from  an  analog  Integrator  and  two  coefficient  poten¬ 
tiometers,  as  shown  in  Figure  12. 

dG 

If  the  integrator  output  is  e  .its  input  is  -  — .  and  thus  the  circuit  solves  the  equation 
o  dt 

de 

-  —  =  +  ore  (80) 

dt  1  o 


FIGURE  11.  RELATIONSHIP  BETWEEN  h(t  -  t  ) 
ANDfjft-tp) 
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The  value  of  K  In  Equation  79  la 

at 


2 

Substituting  -  for  K  In  Equation  79  gives 


K _ : _ - _ ^ 

.00  t 

r  r  f  o  2  2o(t-t 

I  ‘i  I  Qt  e  ( 

J-oo  J -00 


)dt 


(88) 


p((i>  ,  t  )  - -  E  (O),,  t  ) 

O  O  VQt  o  1  O 


(89) 


Figures  13(a)  and  (b)  are  sketches  of  (j(t  -  y  and  jFj(u))|,  the  associated  time  and  frequency  windows 
used  In  the  process. 


|G{w) 


(b) 


FIGURE  13.  ASSCX;iATED  TIME  AND  FREQUENCY  APERTURES,  (a)  Exponential  time 
aperture,  (b)  Frequency  aperture  or  filter  function. 


The  time  aperture  Is  an  exponential  which  gives  greatest  weight  to  the  characteristics  of  f(t)  at  t  -  t  , 

o 

with  the  weight  decreasing  exponentially  as  t  -  t  becomes  larger. 

o 

Substituting  Equations  87  and  68  into  Equation  70  gives 


P(w  .  t  ) 
o  o 


r  t 

2  r-  t 

2I 

/<  0  af(t-t  ) 

fro  ot(t-t  ) 

1  e  °  f(t)  cos  u>  tdt 

+  1  e  °  f(t)  sin  IP  tdt 

-•'-00  °  . 

L  -^-co  °  J 

J 

(90) 


which  Is  the  expression  used  for  evaluating  p((>>  .  t  )  wlien  first  order  low-pass  filters  are  used  for  filters 

o  o 

A  and  B  of  Figure  10. 

Since  the  various  frequency  components  of  seismic  surface  waves  propagate  at  different  velocities,  the 
Initial  part  of  the  received  wave  will  contain  only  those  frequency  components  which  propagate  at  the  higher 
velocities,  and  the  more  slowly  propagating  frequencies  will  arrive  later.  Measurement  of  this  dispersion 
of  the  arrival  times  for  the  different  frequency  components  is  an  important  part  of  analyzing  the  wave.  The 


28 


ImtHut*  of  Scimct  and  Tochnology 


Th«  Univardty  of  Michigan 


exponential  time  aperture  of  Figure  13(a)  ia  well  aulted  (or  thia  type  of  work,  since  it  gives  major  emphasis 
to  the  frequency  characteristics  of  the  portion  of  f(t)  which  is  just  entering  the  window. 

An  expression  resembling  Equation  90  has  been  defined  as  the  e.m.p.  (exponentially  mapped  past)  power 
spectrum  by  Dr.  Joseph  Otterman  |S|.  The  e.m.p.  power  spectrum  is  defined  as  follows; 


r  f(x)e  cos  u)(t  -  x)  dx 

==  2 
+  a 

r  r  f(x)e  alnuJ(t  -  x)dx 

_  J-flO 

L  J-QO 

(91) 


In  order  to  permit  comparison  of  this  expression  with  Equation  90.  we  show  that  cos  uix  may  be  sub¬ 
stituted  for  cos  (U(t  -  X).  and  sin  (Ox  may  be  substituted  for  sin  u)(t  -  x)  In  Equation  91  without  changing  the 
value  of  P^(<0).  Equation  91  can  then  be  written  in  complex  form: 


P^fco) 


"“I. 


f(t)e  cos  a)(t  -  x)dx  +  j  f  f(x)e  sin  a)(t  -  x) 

J  -00 


■“1/ 

(<x)e' 

1  J-ao 

ol  r® 

-  Of  f 

f(x)e' 

1  J  -00 

2I  M 

«  Of  e  1 

i*  OO 

* 

1  ^ 

-00 

-a{t-x) 


( cos  u)(t  -  x)  +  j  sin  u)(t  -  x)l  dx 


-“(t-x)ej‘^(t-x)dx 


-0(1- X)  -juw  . 
e  e  dx 


(92) 


Since  |e  =1,  this  factor  may  be  deleted  from  Equation  92  and  we  can  write 


/>  t 

2  f 

“  -t 

2 

p^(W)  =  or 

f  f(x)e-“<‘-’‘>e-^‘"dx 

J-QO 

-1 

r  f(x)e  cos  dx 

-J-QO 

[/: 


f(x)e  sin  utxdx 


(93) 


Comparison  of  Equation  93  with  Equation  90  reveals  that  since  the  corresponding  integrals  are  Identical, 


P(«  .  t  )  =  —  p  (w) 
0  O  TTO  Of 


(94) 


where  p(a)  ,  t  )  Is  the  time-varying  power  spectrum  defined  by  Equation  90,  and  p  ((O)  is  the  e.m.p.  power 
0  0  0! 

spectrum  defined  by  Equation  91. 
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Also,  according  to  Equation  89 

o  0  irot  0  10 

and  comparison  of  this  with  Equation  94  above  shows  that 

E^  =  P^(U))  (98) 

Thus,  the  output  voltage,  E  ,  of  the  analog-computer  circuit  directly  represents  the  e.m.p.  power  spec- 
o 

trum  as  defined  by  Otterman,  but  this  voltage  must  be  multiplied  by  the  factor  l/tro  to  obtain  the  time- 
varying  power  spectrum  defined  by  Equation  90. 

We  recall  from  the  discussion  following  Equation  30  that  the  derivation  of  p(u)  ,  t  )  is  based  on  a  power 

o  0 

spectrum  which  Includes  negative  as  well  as  positive  frequencies,  and  that 

P(a>  ,  t  )  =  2p(w  .  t  )  (96) 

O  O  0  o 

where  P(<»^,  t^)  is  the  power  density  at  for  time  t^,  and  where  the  frequency  spectrum  of  f(t)  is  considered 
to  exist  for  positive  values  of  <t)  only. 

2.2.3.4.  Change  of  Bandwidth  with  Frequency.  When  the  analog-computer  cicuit  using  first-order  low- 
pass  filters,  as  described  in  Section  2.2.3.3,  is  used  for  measuring  the  power  or  energy  density  at  low  fre¬ 
quencies  such  as  1  cps,  the  filter  passband  must  be  much  narrower  than  1  cps  in  order  to  give  reasonable 
frequency  resolution.  At  higher  frequencies,  such  as  100  cps,  this  extremely  narrow  bandwidth  is  unsatis¬ 
factory  because  a  large  number  of  separate  frequencies  would  have  to  be  investigated  to  cover  a  given  fre¬ 
quency  range.  Consequently,  the  filter  bandwidth  was  changed  for  each  interrogation  frequency.  It  was  ex¬ 
perimentally  determined  that  for  seismic  data,  satisfactory  frequency  resolution  could  be  obtained  by  making 
the  filter  bandwidth  between  the  half-power  points  equal  to  2,5%  of  the  interrogation  frequency  Thus, 
for  the  first-order  low-pass  filters, 


or  solving  for  a. 


bandwidth  =  2a  =  0.025a) 


1 


where  =  a^/2ir 


a  =  0.0125a)j  =  0.0257rF| 


This  value  for  a  may  be  substituted  in  Equation  89  (repeated  here). 


p(a)  ,  t  )  =  — E  (0),,  t  ) 
0  0  xa  0  i  o 


(97) 


(98) 


(89) 
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to  give 


p(w  t  )  - - E  (F  ,  t  )  .  -12-  E  (F,.  t  ) 

0  0  A  AA..  2_  o  I  o  2_  o  i  o 
0.025ir  Fj  IT  Fj 


(99) 


Equation  99  U  valid  only  when  the  bandwidth  2a  la  programmed  to  be  2.5%  of  the  Interrogation  angular 
frequency 

Under  these  conditions,  the  frequencies  to  be  Investigated  were  chosen  0.01  decade  apart,  thus  giving 
200  separate  interrogation  frequencies  in  the  two-decade  range  from  1  cps  to  100  cps.  With  a  0.01 -decade 
separation,  the  relationship  between  two  adjacent  Interrogation  frequencies  is 

logu)  =  logic  +0.01  (100) 

n^l  n 


^n+1 

Solving  this  for - gives 

0) 

n 


(C 

n 


0.01 

10  =  1.0233 


(101) 


Thus,  the  increase  from  any  given  interrogation  frequency  to  the  next  higher  is  2.33%,  and  since  the 
bandwidth  at  each  Interrogation  frequency  is  2,5%,  the  frequency  spectrum  will  be  adequately  covered,  with 
no  gaps  between  adjacent  interrogation  frequencies. 

If  only  positive  frequencies  are  considered  as  making  up  the  seismic  signal,  then  Equation  99  must  be 
multiplied  by  two. 

Combining  Equations  9li  and  99, 


P(‘-o-to)  =  2p(u>„.V~E_,(F,.y 

X  F, 


or,  since  u)  =  w.  and  ic,  =  2irF,,  P  is  also  a  function  of  F,  and  t  ,  and 
o  i  i  1  1  o 


(102) 


TT  F, 


0  10 


(103) 


2.2.3.S.  Computation  of  the  Energy-Density  Spectrum.  If  the  output  of  the  adder  of  Figure  10  is  inte¬ 
grated,  the  result  Is 


"7  =  -  p<“o-  V"* = V"* 


(104) 
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And,  since  the  energy  density  W .  is  given  by 

d 


a  o  o 


/' 


P(w  .  t  )dt 
o  o 


(105) 


(assuming  f(t)  =  0  for  t  <  0) 


(106) 


If  we  substitute  0,02SirF,  for  a  and  designate  E.  and  W,  as  functions  of  the  two  variables,  F,  and  t  , 
1  7  d  1  o 


we  obtain 


(108) 


3 

RISUITS 

Since  the  time-varying  power-density  spectrum  P(F,.  t  ),  discussed  in  Section  2,  is  a  function  of  the 

1  o 

two  variables  F,  and  t  ,  there  are  two  basic  ways  in  which  this  data  could  be  plotted.  Frequency,  F  .  could 
to  1 

be  held  constant  and  the  power  density  P  could  be  plotted  as  a  function  of  time  t^;  or  plots  could  be  made  of 
power  density  versus  frequency  for  various  fixed  values  of  time. 

Figure  14(a)  is  an  illustration  of  the  first  type  of  plot,  applied  to  a  seismometer  signal  produced  by  an 
earthquake.  These  curves,  with  the  exception  of  the  bottom  one,  show  the  product  power  density  and  freq¬ 
uency  at  various  frequencies  as  a  function  of  time.  The  bottom  curve  is  a  plot  of  instantaneous  power  versus 
time  for  the  signal  being  analyzed,  and  is  included  to  serve  as  a  time  reference.  Note  that  although  signals 
were  received  immediately  after  the  sixth  second,  according  to  the  bottom  curve,  the  individual  curves  for 
the  various  frequencies  from  0.5  cps  to  1.41  cps  show  no  significant  activity  before  the  ninth  second.  Fig¬ 
ures  14(b)(c)(d)  and  (e)  are  continuations  of  Figure  14(a)  for  increasingly  higher  interrogation  frequencies. 

As  the  frequency  becomes  higher,  the  starting  time  becomes  earlier  until  it  corresponds  to  the  starting 
time  of  the  seismic  wave  itself. 

Another  interesting  phenomenon  is  the  separation  of  the  trace  into  two  separate  parts,  one  starting  at 
0  seconds  and  the  other  at  about  20  seconds,  for  frequencies  above  7  cps. 

The  second  method  of  presenting  data  (i.e.,  power  density  versus  frequency  for  fixed  values  of  time) 
was  slightly  modified  by  making  the  ordinate  the  product  of  energy  density  and  frequency  instead  of  power 
density. 
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Results  of  this  type  are  shown  in  Figure  IS.  The  curves  labeled  1  through  5  represent  the  signal's 
total  energy-denaity-frequency  product,  as  a  function  of  frequency,  for  times  of  11,  22,  33,  44,  and  55  sec¬ 
onds  from  "zero"  time  shown  on  Figure  14.  Thus,  curve  1  represents  frequency-energy-density  product 
as  a  function  of  frequency,  for  the  first  11  seconds  of  the  signal;  curve  2  is  for  the  first  22  seconds,  and  so 
on.  Finally,  curve  5  shows  the  frequency-energy-density  product  versus  frequency  for  the  entire  seismic 
signal.  The  separation  between  the  curves  at  any  specific  frequency  is  proportional  to  the  average  power 
density  for  the  tinne  Interval  bounded  by  the  two  curves.  A  logarithmic  frequency  scale  was  used  for  Fig¬ 
ure  15  in  order  to  give  uniform  separation  between  the  separate  interrogation  frequencies,  which  are  spaced 
0,01  decade  apart.  (See  Section  2,2.3.4.) 

It  can  be  readily  shown  that  if  the  curves  of  Figure  15  are  extended  to  cover  all  frequencies  from  0  to  <xi, 
the  area  under  the  highest  curve  (number  5)  is  proportional  to  the  total  energy  of  the  seismometer  signal. 
The  area  under  the  curve  is  given  by 


A  = 


(109) 


For  Figure  15, 

y=KjFiWd 


(110) 


where  K  is  a  scaling  constant,  F  is  the  frequency,  and  is  the  energy-density  spectrum  in  watt-seconds 
per  radian  per  second.  Also  for  Figure  15, 


and 


dx  = 


10  F 


(111) 


(112) 


Substituting  into  Equation  109, 


K.K.  -oo 


r 

log  10  J 


WjdF 

d 


(113) 


but  the  total  energy  in  the  signal,  by  the  Fourier  Integral  energy  theorem,  is 


W«f  W.dui«2irf  W.dF 

Jo  **  Jo  '■ 


(114) 
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and,  from  Equation  113, 


Therefore, 


1 


A  log  10 

-4r 


w 


2ir  log  10 
e 


A  =  K  A 
3 


(115) 


(116) 


Thiia,  aside  from  the  fact  that  FW^  is  readily  available  from  the  computer  circuit,  plotting  this  quan¬ 
tity  versus  log  F  instead  of  versus  log  F  has  the  additional  advantage  that  the  area  under  the  curve  is 

proportional  to  total  energy.  (This  would  also  be  true  for  a  plot  of  W .  versus  F,  but  would  not  be  true  for 

d 

W.  versus  log  F.) 
d 

It  is  clear  that  the  area  under  any  one  of  the  curves  of  Figure  15  is  proportional  to  the  total  energy  of 
the  signal  up  to  the  time  represented  1^  that  curve.  The  curves  of  Figures  14  and  15  represent  specific 
results  for  a  particular  seismic  wave.  During  the  investigation,  about  45  separate  tape-recorded  seismic 
records  were  processed  in  the  manner  indicated.  A  detailed  discussion  of  the  seismologlcally  significant 
results  of  this  investigation  is  being  deferred  to  a  later  report,  since  the  present  report  is  concerned  with 
theory  and  techniques  only. 


4 

CONCLUSIONS  AND  RECOMMENDATIONS 


4.1.  CONCLUSIONS 

(a)  For  obtaining  useful  information  from  a  frequency  analysis  of  a  seismic  signal,  the  formal  methods 
which  analyze  the  entire  signal  in  terms  of  constant-amplitude  sinusoids  of  various  frequencies  are  unsatis¬ 
factory.  Instead,  some  method  which  leads  to  a  "time-varying"  spectrum  must  be  used. 

(b)  Any  "time-varying"  spectrum  so  obtained  necessarily  involves  some  approximations,  since  any 
analysis  of  a  signal  in  terms  of  amplitude- modulated  sinusoidal  waveforms  is  not  unique,  but  depends  upon 
the  characteristics  of  the  filter  or  time  window  used. 

(c)  In  obtaining  a  time-varying  power  spectrum  by  any  technique,  the  time  resolution  and  frequency 
resolution  are  interrelated  in  such  a  manner  that  an  improvement  in  one  results  in  a  degradation  of  the 
other.  They  can  both  be  improved  somewhat  by  optimum  choice  of  time  and  frequency  apertures,  but  the 
basic  problem  arises  from  the  Interrelated  nature  of  the  time  and  frequency  domains,  and  not  from  equip¬ 
ment  limitations. 
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(d)  Three  possible  techniques  for  obtaining  s  time-varying  power  spectrum  are: 

1,  Weighted  Fourier  analysis,  using  a  moving  "time  window"  as  the  weighting  fonction. 

2,  Use  of  variable  filters  or  "frequency  apertures"  and  measuring  the  amount  of  power,  as  a  func¬ 
tion  of  time,  which  passes  through  the  filter  for  various  frequencies. 

3,  Translation  of  the  frequency  spectrum  of  the  seismic  wave  by  heterodyning  techniques  in  such 
a  manner  that  each  portion  of  the  spectrum  in  turn  may  be  passed  through  a  fixed-frequency 
filter. 

Theoretically,  these  three  techniques  are  equivalent,  and  no  one  method  can  give  results  which  are  superior 
to  those  obtained  from  the  others.  Choice  of  the  technique  to  be  used  should  be  based  on  available  equipment 
and  its  limitations. 

(e)  Although  it  offers  no  theoretical  advantage  over  methods  using  filters,  the  analog  conq>utlng  tech¬ 
nique  described  in  this  report  la  by  its  flexibility  probably  is  the  most  satisfactory  for  developing  and  testing 
various  methods.  After  an  optimum  method  of  analyzing  seismic  data  has  been  found,  it  would  probably  be 
advantageous  to  build  special  apparatus  for  routine  analysis  rather  than  continue  using  an  analog  computer. 
Such  apparatus  would  probably  be  less  expensive  to  operate  and  would  perform  the  analysis  much  more 
rapidly  than  the  analog  computer  setup, 

(f)  Seismic  spectra  in  general  are  not  smooth,  but  contain  a  number  of  definite  peaks.  In  order  to  re¬ 
solve  these  peaks,  a  filter  bandpass  of  about  2  1/2%  of  the  center  frequency  appears  necessary.  No  signif¬ 
icant  improvement  in  resolution  is  obtained  by  using  narrower  filters. 

4.2.  RECOMMENDATIONS 

Recommendations  for  future  work  may  be  divided  into  two  categories:  (a)  optimization  of  frequency- 
analysis  techniques;  and  (b)  development  of  other  methods  of  analyzing  seismic  data.  The  principal  object 
of  future  investigation  would  be  to  develop  methods  of  analyzing  seismic  waves  which  lead  to  accurate  deter¬ 
mination  of  the  nature  and  location  of  the  sources  and  the  physical  characteristics  of  the  transmission  path. 
Investigation  of  the  spectral  characteristics  of  a  seismic  wave  is  only  one  of  several  reseach  programs 
which  may  lead  to  significant  results. 

4.2.1.  FURTHER  DEVELOPMENT  AND  REFINEMENT  OF  FREQUENCY-ANALYSIS  TECHNIQUES.  We 
recommend  that  further  research  be  performed  on  the  Items  below. 

(a)  Theoretical  and  experimental  investigations  to  find  a  technique  which  gives  the  best  compromise 
between  time  and  frequency  resolution  when  applied  to  seismic  data. 

(Although  the  time  and  frequency  apertures  may  not  be  independently  chosen  since  the  choice  of  one 
automatically  specifies  the  other,  it  seems  possible  to  find  a  system  which  gives  optimum  results  in  terms 
of  time  resolution  and  frequency  resolution  considered  simultaneously.  A  particular  frequency  aperture  to 
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be  tested  is  one  with  a  shape  corresponding  to  the  bell-shaped  Gaussian  distribution  curve.  Since  the  Fourier 
transform  of  this  curve  is  also  a  Gaussian  distribution  curve,  the  time  aperture  in  this  case  would  be  iden¬ 
tical  in  shape  with  the  frequency  aperture.) 

(b)  Attempting  to  find  a  graphical  presentation  type  of  a  time-varying  spectrum  which  will  result  in  a 
"signature"  characteristic  of  the  type  of  seismic  wave  being  analyzed. 

(c)  Investigating  the  possibility  of  accurately  measuring  the  arrival  times  of  the  principal  frequency 
components  of  the  seismic  wave  by  first  using  a  narrow  filter  to  determine  the  center  frequencies  of  peaks 
in  the  spectrum,  and  then  using  wideband  filters  to  measure  the  arrival  time  of  these  frequencies. 

(d)  Investigating  the  feasibility  of  playing  the  tape-recorded  seismic  signal  both  forward  and  backward  to 
measure  more  accurately  the  arrival  time  of  each  frequency  component. 

4.2.2.  DEVELOPMENT  OF  METHODS  OTHER  THAN  FREQUENCY  ANALYSIS  FOR  ANALYZING  SEIS¬ 
MIC  DATA.  Other  methods  of  analysis  which  could  be  investigated  are  the  method  of  time-varying  autocorre¬ 
lation  functions  and  that  of  time-varying  amplitude-distribution  functions. 

(a)  The  autocorrelation  function  of  a  function  of  time  is  useful  not  only  in  power-spectra  determination, 
but  also  in  detecting  periodic  signals.  Any  periodic  signals  which  were  concealed  In  seismic  data  would 
probably  appear  only  during  certain  time  intervals.  Thus,  techniques  for  determining  a  time-varying  auto¬ 
correlation  function  could  be  developed  and  applied  to  seismic  data. 

(b)  A  signal  which  has  uniform  statistical  characteristics  for  any  interval  of  time  (stationary  function) 
is  often  described  in  terms  of  its  amplitude-distribution  function.  Since  the  statistical  characteristics  of  a 
seismic  wave  are  not  uniform  in  time,  the  amplitude-distribution  function  as  measured  for  one  time  interval 
may  be  considerably  different  from  that  determined  for  the  preceding  or  following  time  interval.  Analog 
computing  techniques  similar  to  those  used  for  time-varying  power-spectra  measurement  could  be  used  to 
determine  a  time-varying  amplitude-distribution  function.  Whether  or  not  such  a  function  would  be  useful 
for  describing  seismic  signals  in  terms  of  their  source  and  transmission  path  would  be  determined  experi¬ 
mentally. 
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Appondix  A 

SEISMIC- SION AL  CORRECTION 

In  the  following,  the  notation  f(t)  designates  a  signal  produced  by  a  seismometer  with  a  flat  response, 
and  f'(t)  represents  the  output  of  a  seismometer  which  is  sensing  the  signal.  In  general,  f'(t)  differs  from 
f(t)  in  the  seismometer's  response  characteristics.  The  seismometer  output,  f'(t),  had  been  recorded  on 
magnetic  tape,  using  an  FM  recording  system,  and  the  analysis  was  performed  on  the  signal  played  back 
from  the  tape.  No  significant  frequency  distortion  resulted  from  the  tape-recording  and  reproducing  proc¬ 
esses,  but  a  slight  difference  between  the  center  frequency  of  the  playback  demodulator  and  the  carrier 
frequency  of  the  recorded  signal  generated  a  d-c  bias  voltage  which  was  not  originally  present  at  the  seis¬ 
mometer  output.  In  order  to  obtain  an  accurate  reproduction  of  f(t)  for  frequency  analysis,  both  the  d-c 
introduced  by  the  record-reproduce  process,  and  the  spectral  distortion  introduced  by  the  seismometer  re¬ 
sponse  were  removed. 

If  the  complex  frequency  spectrum  of  f(t)  is  F(s),  and  the  seismometer  response  is  G(s),  the  output  of 
the  seismometer  is  given  by 

F-(s)  =  F(s)G(s)  (117) 


Since  the  record- reproduce  process  adds  a  d-c  bias,  the  output  of  the  playback  demodulator  is  given  by 

F"(s)  =  F'(s)  +  bias  =  F(s)C(s)  +  bias  (118) 

Equation  118  can  be  solved  for  F(s)  to  give 


F(s)  = 


F"(8)  -  bias 
G(8) 


(119) 


and  the  time  response  corresponding  to  F(s)  then  represents  f(t). 

The  recording  process  is  illustrated  by  the  upper  block  diagram  of  Figure  16,  The  lower  block  diagram 
of  this  figure  illustrates  the  arrangements  of  the  correcting  filters  in  the  playback  setup.  The  FM  signal 
from  the  tape-playback  unit  enters  the  demodulator,  where  it  is  converted  to  a  signal  corresponding  to 
F'(s)  or  f'(t).  The  filter  block  with  the  transfer  characteristic  ^  removes  any  d-c  bias  from  the  output 
of  the  demodulator,  and  the  block  labeled  1/G(s)  removes  the  effect  of  the  seismometer  response  from 
the  signal.  The  output  of  the  Inverse  seismometer  circuit  is  then  approximately  f(t),  the  error  depending 
upon  how  accurately  the  Inverse  seismometer  response  is  represented  by  the  1/G(8)  block. 

Average  frequency-response  characteristics  of  two  types  of  seismometers  (Benioff  and  Willmore)  are 
shown  in  Figure  17.  These  curves  represent  the  average  of  a  number  of  curves  obtained  from  shaketable 
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teiU.  For  these  tests,  the  Input  signal,  as  obtained  from  the  shaketable,  was  a  sinusoid  having  an  rms 
velocity  of  0.01  Inches  per  second.  The  curves  show  the  resulting  rms  voltage  output  of  the  seismometers 
as  a  function  of  frequency.  In  Figure  17,  the  peak  response  of  the  Benloff  curve  Is  about  0.8  volt  at  2.2  cps. 

If  the  Inverse  seismometer  response  fitter  had  a  gain  of  unity  at  this  frequency,  the  voltage  scaling  from 
the  seismometer-filter  combination  would  be  O.g  volt  per  Inch  per  second.  If  the  correcting  filter  has  a 
gain  of  1/0.8  at  this  point,  however,  the  output  scaling  of  the  seismometer-corrector  combination  Is  1  volt 
per  Inch  per  second  over  the  frequency  range  for  which  the  corrections  are  valid. 

A.I.  BENIOFF-SEISMOMETER  RESPONSE  CORRECTION 

The  response  curve  for  the  Benloff  seismometer  Is  redrawn  In  Figure  18.  Below  It  a  frequency- response 
curve  Is  plotted  as  measured  from  an  analog-computer  mechanization  which  was  used  to  approximate  the 
Inverse  seismometer  response. 

The  resultant  frequency  response  of  the  seismometer-filter  combination  Is  the  product  of  the  two 
curves;  It  Is  also  plotted  In  Figure  18. 


As  shown,  the  resultant  droops  at  both  the  high-  and  low-frequency  ends,  but  since  the  range  of  the  In¬ 
vestigation  was  from  O.S  to  SO  cps,  with  the  major  portion  of  the  return  well  within  these  limits,  it  was  felt 
that  this  droop  was  permissible. 


FREQUENCY  (cps) 


FIGURE  18.  BENIOFF-SEISMOMETER  CORRECTION  CURVES 
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The  transfer  function  G(s)  between  the  velocity  of  the  shaketable  and  the  output  of  the  Benloff  seismom¬ 
eter  Is  approxlnoated  by 


8*  Zt 

1 

2  (D 
b)  n 
n 


where  K  =  a  gain  constant 


=  the  natural  angular  frequency 


S  =  a  damping  factor 


The  inverse  of  this  function  Is 


1  _1 
G(s)  °  K 


2  b)  S 
u>  n 
n 


In  order  to  avoid  difficulties  arising  from  differentiating  high-frequency  noise  and  from  Integrating  low- 
frequency  drift,  the  transfer  function  for  the  Inverse  response  was  changed  to 


_L*i 

G(s)  K 


2„  ,  b>  T,s  +  1 

n  1 


The  mechanization  of  - — ^  is  an  approximate  integrator  which  is  good  for  frequencies  much  greater  than 
’'i®  V 

1/t,  radians  per  second,  and  the  mechanization  of  - Is  an  approximate  differentiator  which  is  valid  for 

1  1  +  T_8 

2 

frequencies  much  less  than  I/t^  radians  per  second.  Rearranging  terms  of  Equation  122  in  order  of  ascend¬ 
ing  "break-point"  frequencies  of  the  approximate  integrator  and  the  approximate  differentiator  gives 


-L-.i  ^—3  . _ ? _ 

C(s)  K  TjS.l  b,^ 

—  n  2  _ 


Figures  19(a)  and  19(b)  show  the  response  curves  that  would  result  if  perfect  Integrators  and  differen¬ 
tiators  were  used  to  generate  the  Inverse  Benloff  response  characteristics  called  for  in  Equation  121.  Here 
the  response  curve  would  start  with  a  gain  of  infinity  for  d-c  (curve  a)  and  follow  a  negative  integrator  slope 
down  with  increasing  frequency  until  it  met  the  added  positive  slope  (curve  b)  of  the  differentiator,  where  the 
gain  would  rise  again  back  to  infinity  for  infinite  frequency.  The  constant  term  of  Equation  121  adds  to  the 
gain  value  and,  at  the  point  of  intersection  of  the  negative  and  positive  gain  slopes,  adjusts  the  value  and 
shape  of  the  minimum  gain  area.  Figures  19(c)  and  19(d)  show  the  response  curves  of  an  approximate  inte¬ 
grator  and  differentiator  which  are  specified  by  Equation  123.  From  inspection  of  Figures  19(c)  and  19(d), 
it  is  apparent  that  the  break-point  frequencies  must  be  moved  up  the  slopes  as  far  as  possible  so  that  the 
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FIGURE  19.  RESPONSE  CURVES  OF  ACTUAL  AND  APPROXIMATE  INTEGRATORS 
AND  DIFFERENTIATORS,  (a)  Integrator,  (b)  Differentiator,  (c)  Approximate  integrator, 
(d)  Approximate  differentiator. 


actual  curve  closely  matches  the  ideal  curve  within  the  frequency  range  of  interest.  This  was  done  and  the 
inverse  Benioff  response  curve  was  followed  almost  perfectly  for  all  teat  frequencies  used;  but  for  the 
actual  runs  the  high-frequency  gain  was  too  much  for  the  electronic  multiplier.  If  the  frequency  of  the  signal 
and/or  noise  is  much  above  100  cps,  the  amplitude  must  be  greatly  reduoxl  in  order  to  keep  the  electronic 
multiplier  from  overloading. 
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High-frequency  response  was  suniressed  enough  to  satisfy  the  requirements  of  the  electronic  multiplier 
by  adding  a  small  feedback  capacitor  around  the  output  amplifier  of  the  Inverse  filter,  as  shown  In  Figure  20. 

In  order  to  shorten  the  time  constant  of  the  d-c  bias  eliminator  (an  approximate  differentiator),  Ita  break 
frequency  was  moved  up,  which  made  the  addition  of  another  approximate  integrator  necessary  to  increase 
the  negative  slope  at  the  low-frequency  end  of  the  Inverse  Benioff  filter.  The  completed  filter  had  a  frequency- 
response  curve  which  accurately  approximated  the  inverse  Benioff  response  characteristics  within  the  fre¬ 
quency  range  of  interest,  and  which  dropped  att  rapidly  for  frequencies  outside  this  range. 

Figure  20  shows  an  analog-computer  diagram  for  the  solution  of  the  inverse  Benioff  response  charac¬ 
teristic.  The  parts  which  were  added  to  compensate  for  the  analog  computer  limitations  are  drawn  with 
dashed  lines. 


FIGURE  20.  ANALOG-COMPUTER  CIRCUIT  FOR  INVERSE  BENIOFF-SEBMOMETER  RESPONSE 


A.2.  WILLMORE-SEISMOMETER  RESPONSE  CORRECTION 

The  Willmore-seismometer  response  and  the  inverse  response  function  used  are  shown  in  Figure  21. 
Since  the  gain  at  the  high-frequency  end  is  only  about  3.3  for  a  maximum,  no  trouble  was  encountered  in  the 
solution. 
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FREQUENCY  (cps) 


FIGURE  21.  WILLMORE-SEISMOMETER  CORRECTION  CURVES 


The  transfer  operator  for  the  Willmore  seismometer  was  approximated  by 


G  (s) 
w 


8  2£ 

- +  —  s  +  1 

2  b) 

W  n 

n 


The  inverse  of  this  operator  is 


(124) 


1 

G  (8) 
w 


(126) 
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(126) 


(127) 


FIGURE  22.  ANALOG-COMPUTER  CIRCUIT  FOR  THE  INVERSE  WILLMORE- 
SEI8MOMETER  RESPONSE  FILTER 


Appendix  B 

ANALOG -COMPUTER  IMPLEMENTATION 

This  appendix  contains  a  more  detailed  description  of  the  analog  computing  circuits  than  that  presented 
in  Section  2.2,3.  of  this  report.  The  main  computing  circuits  were  assembled  on  a  "PACE"  computer  made 
by  Electronic  Associates,  Inc.  The  automatic  programming  was  performed  by  additional  conventional  analog 
equipment  built  in  the  1ST  Analog  Computer  Laboratory.  The  interrogation  frequency  for  each  run  was 
selected  by  a  dlgital-to-analog  converter,  operating  from  a  punched -paper-tape  reader. 
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B.1  ANALOG-COMPUTER  ARRANGEMENT 

Figure  23  is  a  block  diagram  of  the  complete  analog-computer  aetup  used  for  the  frequency  analyala  of 
selamlc  stgnala.  It  conalats  of  the  baalc  computing  block  diagram  of  Figure  10  and  the  blocks  used  for  cor¬ 
recting  the  recorded  seismic  signal,  controlling  the  computer,  and  plotting  the  final  results.  Input  informa¬ 
tion  is  received  from  two  tape-recorder  channels,  shown  at  the  upper  left  in  Figure  23.  Channel  1  is  the  con¬ 
trol  channel  and  channel  2  is  the  signal  or  information  channel. 

The  signal  from  the  demodulator  for  channel  2,  the  Information  channel,  is  r'(t).  As  explained  in  A^ien- 
dix  A,  this  signal  consists  of  f(t)  multiplied  by  the  seismometer  response  and  containing  a  small  d-c  bias. 

The  bias  eliminator  is  a  high-pass  filter  which  removes  any  d-c  component  from  the  signal.  Its  output  is 
r(t),  the  seismic  signal  f(t)  multiplied  by  the  seismometer  response  characteristics.  Switches  SWl  and  SW2 
may  be  positioned  either  to  pass  f'(t)  directly  into  the  computer  circuits,  or  to  insert  the  proper  inverse 
seismometer  filter  and  thus  produce  f(t),  corresponding  to  the  output  of  a  seismometer  with  a  flat  response. 

The  corrected  signal  f(t)  then  enters  the  basic  computing  circuit,  whose  operation  has  been  explained 
qualitatively  In  Section  2.2.3.  In  order  to  explain  the  operating  sequence  and  at  the  same  time  quantify  the 
results  somewhat,  we  assume  that  the  signal  f(t)  Is  replaced  by  a  sinusoid  A  sin  out.  Then,  the  block  marked 
''Multiplier  A"  in  the  upper  channel  multiplies  the  signal  A  sin  out  by  cos  ou^t  from  the  oscillator  block  to  pro¬ 
duce  at  its  output 

A  sin  out  cos  0U|t 
or 

E^  -  y  (sin  (<u  -  0Uj)t  +  sin  (ou  +  Wj)t] 

and  multiplier  B  in  the  lower  channel  multiplies  the  signal  A  sin  out  by  sin  ou^t  from  the  oscillator  to  produce 
at  its  output 

Ej  =  J  [cos  <<*»  -  0Uj)t  +  cos  (OU  +  OUjlt) 

Filter  A  in  the  oipper  channel  and  filter  B  in  the  lower  channel  are  both  low-pass  filters  and,  as  such,  attenu¬ 
ate  all  signals  except  those  of  zero  frequency.  Consequently,  the  sum-frequency  terms  ( sin  (ou  +  ou^ft  and 
cos  (ou  +  0U|)t]  in  the  above  equations  are  attenuated  to  such  a  degree  that  their  amplitudes  may  be  considered 
negligible  at  the  outputs  of  the  filters.  Thus,  the  output  of  filter  A  of  the  upper  channel  is 

A 

Eg  -  j  sin  (ou  -  oUjlt  G(tu  -  oUj) 
and  the  output  of  filter  B  of  the  lower  channel  is 

■  '*’i) 

4  2  I  1 
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The  computer  control  block  operates  together  with  the  digital -to-analog  (D/A)  converter  to  generate  an 
analog  voltage  corresponding  to  a  particular  interrogation  frequency  for  each  computer  run.  This  voltage 
operates  servo-positioned  potentiometers  in  the  oscillator,  in  order  to  set  it  at  the  corresponding  frequency, 
and  in  filters  A  and  B  so  that  the  bandwidth  of  the  filters  will  be  proportional  to  the  oscillator  frequency. 
Multipliers  C  and  D  are  servo  multipliers  used  to  square  the  outputs  of  filters  A  and  B,  respectively.  By 
use  of  a  gain  of  2  at  the  multiplier  Inputs,  the  1/2  is  eliminated  from  the  equation  of  the  output  of  the  filters; 
therefore  the  output  of  multiplier  C  becomes 

Eg  -  A^  sln^  (a>  -  g2(w  -  Wj) 

for  the  upper  channel  and 

2  2  2 
E^,  »  A  cos  Ifi)  -  iii|)t  G  (w  -  oij) 

for  the  lower  channel. 

The  summing  amplifier  adds  these  two  multiplier  outputs  (E.)  and  (E„)  to  produce  at  its  output 

5  6 

E^  =  Eg  +  Eg  =  A*G*(W  -  W|) 

This  result  is  integrated  to  produce  at  the  integrator  output 

2  2 

E  =  j  A  G  (W  -  u)  )dt 

’  Jo  ' 

where  T  is  the  time  of  the  seismic  signal  on  the  magnetic  tape,  or  more  precisely,  the  time  of  the  "On"  sig¬ 
nal  in  the  control  channel. 

The  oscillating  term  has  disappeared  In  the  equation  for  the  output  of  the  adder: 

2  2 

E  =  A  G  (O)  -  ii>,) 
o  1 

However,  it  should  be  pointed  out  that  the  difference  frequency  (b>  -  0)^)  determines  the  filter  response 
G(a)  -  u)^),  and  thus  determines  the  attenuation  of  A,  whereas  interrogation  frequency  (u)|)  controls  the  band¬ 
width  of  the  filter. 

The  control  channel  in  Figure  23  is  operated  by  a  positive  or  negative  voltage  from  demodulator  1.  The 
control  signals  were  recorded  on  a  separate  channel  of  the  magnetic  tape  which  contained  the  seismic  signal 
to  be  analyzed.  A  positive  voltage,  or  "On"  signal,  was  recorded  for  the  duration  of  the  seismic  signal,  A 
negative  voltage,  or  "Ofr'  signal,  was  then  recorded  for  about  fifteen  seconds  after  terminating  the  seismic 
signal.  This  magnetic  tape  was  spliced  into  an  endless  loop  which  could  be  played  continuously. 
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The  operating  sequence  of  the  computer  control  followe.  During  the  playback  of  the  aeiimic  signal,  the 
control'channel  voltage  keeps  the  computer  on  "Operate,"  and  the  seismic  signal  is  investigated  for  frequency 
components  in  the  neighborhood  of  the  interrogation  frequency  to  which  the  oscillator  is  set.  At  the  end 
of  the  seismic  signal,  the  "On"  voltage  from  the  control  channel  changes  to  an  "Off  voltage;  this  change  is 
sensed  by  the  computer  control  and  causes  it  to  operate  a  relay  which  in  turn  switches  the  computer  to  the 
"Hold"  mode.  This  stops  all  Integrators  In  the  computing  circuit  and  maintains  all  analog  voltages  at  their 
values  at  the  end  of  the  run.  The  "Off  voltage  in  the  control  channel  also  causes  the  controller  to  start  a 
timing  integrator  whose  output  Is  connected  to  various  differential  relays,  so  arranged  and  connected  that 
the  following  sequence  is  followed. 

After  a  time  Interval  long  enough  for  the  servo-driven  multipliers,  resolvers,  and  plotter  to  come  to 
rest,  a  relay  operates,  montentarily  dropping  the  pen  on  the  plotter,  and  thereby  plotting  the  value  of  E,^ 
from  the  integrator.  Next,  the  hold  condition  is  momentarily  removed  from  the  digital-to-analog  converter, 
and  the  paper  tape  advances,  reading  the  next  set  of  values  for  the  oscillator  and  the  filters,  and  f^^^ 
for  the  plotting  table).  Next,  the  computer  is  placed  in  the  "Reset"  mode,  which  establishes  all  initial  volt¬ 
ages  required  for  the  problem  at  time  "zero".  The  computer  is  now  ready  to  receive  an  "On"  signal  to  start 
the  next  computation  at  the  beginning  of  the  next  playback  cycle  of  the  tape  loop. 

The  D/A  converter  reads  from  a  punched-paper  tape  on  command  from  the  computer  control,  and  feeds 
this  information  back  to  the  computer  control  in  the  form  of  a  d-c  voltage.  The  plotter  is  an  X-Y  plotting 
table  having  an  external  pen-drop  circuit.  The  X  is  positioned  by  a  voltage  proportional  to  the  interrogation 
frequency  (w^);  the  Y  is  positioned  by  the  value  of  the  voltage  at  the  output  of  the  integrator  E,^  at  the  end  of 
each  integration.  Thus  a  line  drawn  through  these  integration  end  points  produces  a  curve  showing  energy 
versus  frequency  for  the  entire  signal.  More  detailed  descriptions  of  the  individual  blocks  of  Figure  23 
follow. 

B.1.1.  THE  OSCILLATOR.  The  oscillator,  which  generates  sin  u^^t  and  cos  w^t,  is  basically  a  three- 
amplifier  oscillator  loop  where  two  of  the  amplifiers  are  Integrators.  This  arrangement  is  shown  schemat¬ 
ically  in  Figure  24. 

B.1.2.  THE  ELECTRONIC  MULTIPLIER.  Multipliers  A  and  B  of  Figure  23  were  actually  two  channels 
of  a  sln^e  electronic  multiplier.  This  is  a  four-quadrant  multiplier  in  that  any  combination  of  mathematical 
signs  of  input  variables  will  always  produce  the  correct  sign  on  the  output.  The  bandwidth  of  this  multiplier 
is  over  100  cps. 

B.1.3.  THE  LOW- PASS  FILTER.  Filters  A  and  B  of  Figure  23  were  simple  first-order  low-pass  units 
assembled  from  analog  components,  as  shown  in  Figure  25. 
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FIGURE  25.  ANALOG-COMPUTER  SCHEMATIC  OF  A  LOW-PASS 
FILTER 


If  the  integrator  output  is  E  ,  its  input  is  -E  ,  and  the  circuit  solves  the  equation 
o  o 

E  +  otE  »  -otE. 
o  o  i 

Obtaining  the  Laplace  transform  of  both  sides  of  Equation  80  and  assuming  that  E  (0)  is  zero, 

0 

sE  (8)  oE  (8)  ■  'aE^(8) 

0  0  1 


Eo<S) 

Ej(e) 


cr 

8 't'  or 


1 

1  +  BT 


(80) 


(128) 


(129) 
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where  t  ■  1/a. 

If  the  filter  circuit  shown  in  Figure  2S  is  to  have  unity  d-c  gain,  it  is  inq>ortant  that  the  potentiometer 
labeled  (a)  in  the  feedback  path  of  the  integrator  have  exactly  the  same  setting  as  the  (a)  potentiometer  in  the 
input  circuit.  As  with  the  previously  described  oscillator,  these  potentiometers  are  set  by  a  servo  multiplier 
so  that  the  value  of  a  la  a  function  of  the  interrogation  frequency  U)^/2ir.  Any  mismatch  between  these  two 
potentiometers  will  affect  the  d-c  gain  of  the  filter  and  cause  ah  apparent  error  in  the  filter  response. 

Figure  26  shows  a  rearrangement  of  the  filter  so  that  and  are  summed  before  multiplying  by  a 
instead  of  multiplying  each  by  a  and  then  summing  at  the  input  of  the  integrator.  This  arrangement  insures 
a  gain  of  unity. 


FIGURE  26.  ANALOG-COMPUTER  SCHEMATIC  OF  A  LOW-PASS  FILTER  USING 
ONLY  ONE  (a)  POTENTIOMETER 


B.1.4.  THE  SQUARING  MULTIPLIER.  Multipliers  C  and  D  In  Figure  23  were  high-speed  servo  units 
used  for  squaring  the  signals  from  filters  A  and  B,  respectively,  before  they  were  added.  Althou^  electronic 
multipliers  or  function  generators  could  have  been  used  for  this  purpose,  the  servo  multipliers  were  used 
because  they  were  readily  available  and  their  speed  was  compatible  with  the  frequency  spectra  of  the  slgnala 
coming  from  the  low-pass  filters. 

A  diagram  of  the  squaring  circuit  is  shown  in  Figure  27.  The  output  of  this  circuit  is  always  positive, 
regardless  of  the  polarity  of  the  input. 

One  particular  advantage  of  the  use  of  a  servo  multiplier  for  performing  the  squaring  function  is  the 
lack  of  any  residual  output  voltage  when  the  input  is  zero.  When  diode  function  generators  were  used  for 
squaring,  considerable  difficulty  was  encountered  in  reducing  the  drift  resulting  from  output  voltages  pro¬ 
duced  by  the  function  generators  at  zero  input. 

B.1.5.  THE  COMPUTER  CONTROL  AND  DIGITAL-TO-ANALOG  (D/A)  CONVERTER.  The  computer  con¬ 
trol  circuits  are  shown  in  Figure  28.  When  an  off  signal  or  minus  voltage  is  received  from  the  control  channel ,  a 
switching  amplifier  (labeled  "H.G.”  in  Figure  28)  switches  the  computer  to  the  "Hold"  mode  of  operation. 
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FIGURE  27.  SERVO  MULTIPLIER  USED  FOR  SQUARING 


thus  retaining  all  voltages  present  at  this  point  of  the  computation.  The  operation  of  a  switching  amplifier 
such  as  this  is  based  on  the  action  of  a  hi^-gain  amplifier  with  two  different  diode-feedback  circuits.  When 
the  input  voltage  E^  is  positive,  the  output  voltage  E^  tends  to  be  negative,  but  the  diode  Lj  produces  a  feed¬ 
back  path  of  approximately  zero  Impedance;  hence  *  0  and  sensitive  relay  SK^  is  unenerglzed  as  shown. 
When  Ej  is  positive,  E^  tends  to  go  positive  and  diode  L^  does  not  conduct.  Under  these  circumstances,  the 
feedback  path  must  be  through  diode  L^.  This  feedback  path  is  esUbllshed  for  a  positive  E^  when  th^voluge 
at  Junction  E^  (between  R^  and  R^)  is  approximately  zero.  This  action  occurs  only  when  ® 

Rj  2 

Thus,  the  output  of  the  switching  amplifier  is  either  0  or  +100  —  volts,  depending  upon  the  polarity  of  the 

**2 


input  voltage.  With  the  proper  ratio  of  R^  to  R^,  the  minus  "Off  voltage  produces  a  plus  E^  of  sufficient 
amplitude  to  energize  relay  SKj  and  thereby  switches  the  computer  to  the  "Hold"  mode.  The  RC-tlme-con- 
stant  filter  placed  between  the  contacts  of  SK^  and  the  hold  relay  reduces  the  off  signal  to  a  pulse  and  thus 
releases  the  hold  command  so  that  the  computer  can  cycle. 


The  hold-command  pulse,  coming  through  network  RC  from  the  switching  amplifier,  also  energizes  relay 
Kj,  and  the  "feedback"  voltage  from  the  computer  hold  circuits  keeps  energized.  Relay  switches  volt¬ 
age  t  into  integrator  1  and  removes  the  short-circuit  from  the  integrator's  feedback  path,  thus  allowing  the 
integration  of  t.  The  voltage  output  (t)  of  integrator  1  rises,  thereby  supplying  relays  SK^  and  SK^  with  an 
increasing  voltage.  The  most  sensitive  relay  was  chosen  for  SK^;  thus  it  operates  first,  produces  a  "pen- 
drop"  signal  through  the  RC  network,  and  causes  the  X-Y  plotter  to  produce  a  dot.  Relay  SK^,  operates 
next;  it  removes  the  "Hold”  voltage  from  the  D/A  converter;  allows  it  to  read  the  next  b)^;  and  also  puts  the 
computer  on  the  "Reset"  mode,  thus  removing  the  computer  Hold  and  de-energizing  relay  K^.  When  the  com¬ 
puter  goes  on  "Reset,"  the  Integrator  output  t  goes  to  zero,  and,  de-energizing  SK^  and  SK^,  re-establishes 
the  D/A  hold  and  sets  the  computer  for  the  "Operate"  command. 
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